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Abstract 

I develop a renormalization-group blocking framework for lattice QCD with staggered fermions. 
Under plausible, and testable, assumptions, I then argue that the fourth-root recipe used in numer- 
ical simulations is valid in the continuum limit. The taste-symmetry violating terms, which give 
rise to non-local effects in the fourth-root theory when the lattice spacing is non-zero, vanish in 
the continuum limit. A key role is played by reweighted theories that are local and renormalizable 
on the one hand, and that approximate the fourth-root theory better and better as the continuum 
limit is approached on the other hand. 
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I. INTRODUCTION 



Lattice QCD simulations with staggered fermions [l|, 0, [s], 0, SJS] have been pro- 
ducing remarkably accurate predictions of various hadronic observables [9j. The staggered- 
fermion field has only one component per color per lattice site, making the numerical com- 
putations relatively cheap, as well as a non-anomalous U{1) chiral symmetry in the massless 
limit, which is important for the phenomenology of the light-quark sector. 

All staggered-fermion simulations with three fiavors of light quarks make use of the fourth- 



root recipe [10[- The up, down, and strange quarks are each represented by a staggered field 
with a different bare mass. But normally a single staggered field generates four quark species, 
or "tastes," in the continuum limit. The four tastes do have equal renormalized masses 
thanks to the lattice staggered-fermion symmetries [S].^ In order to remove the excessive 
degrees of freedom one takes the fourth root of the staggered-fermion determinant.^ The 
fourth-root recipe defines a renormalizable theory which, to all orders in perturbation theory, 
reproduces a local, unitary theory with the correct number of light quarks in the limit of a 
vanishing lattice spacing. (It is assumed [l2| that the staggered theory without the fourth 



root behaves as expected in perturbation theory, which is very plausible [ll|, |l3|.) Non- 



perturbatively, the validity of the fourth-root recipe is a non-trivial issue which has been 



the subject of much debate [IJ]. For a recent review, see Ref. [11 . 

In a formal expansion in the lattice spacing, the massless staggered action splits into 
marginal terms that have a U{4) taste symmetry and irrelevant terms that break the sym- 
metry explicitly. Because of the absence of an exact four-fold degeneracy in the spectrum 
of the staggered Dirac operator, the fourth-root theory is non-local, and not unitary, at 



any non-zero lattice spacing [15|. The question that must be addressed is whether the 
non- locality disappears, and unitarity is recovered, in the continuum limit. 

First, a degree of control over the infra-red behavior must be maintained, and I will 



assume that the quark masses are all positive [16|, ll7|, Il8|, Il9| . In a nut-shell, the following 
tentative reasoning summarizes how the fourth-root recipe might be valid in the continuum 
limit: The taste- violating effects of staggered fermions arise from irrelevant operators. These 
effects should scale to zero in the continuum limit like a positive power of the lattice spacing. 
Hence exact four-fold taste degeneracy is recovered in the continuum limit. The effective 
low-energy Dirac operator attains the form D ® 1, where D is a local operator that carries 
no taste index, and where 1 is the four-by-four identity matrix in taste space. The fourth 
root of det{D 1) is det{D). This fourth root is analytic, leading to locality and unitarity.^ 
Before making any concrete claims, a framework is needed where statements about 
the continuum limit will be well defined to begin with. The natural tool for the task is 
renormalization-group (RG) block transformations, which have already been used in the 
free theory j20|. A coarse-lattice theory is obtained via n blocking steps, starting from a 
fine-lattice theory which, in the case at hand, contains staggered fermions. The process is re- 
peated with more and more blocking steps, but the coarse-lattice spacing (in physical units) 



^ In principle, one could account for the up, down, strange, and charm quarks by the four tastes of a single 
staggered field with a general staggered mass term _8]. In practice, this is not done. See Ref. fill for a 
discussion of the reasons. 

^ In the isospin limit, the up-down sector is represented by a square root of a staggered determinant with 

the common light quark mass. 
^ The analytic continuation to Minkowski space must be performed after the continuum limit. 
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is held fixed. With each additional blocking step the initial fine-lattice spacing gets smaller, 
while the bare parameters are adjusted to maintain constant physics. In the limit n — >■ oo, 
one obtains a (well-defined!) coarse-lattice theory that generates a set of continuum-limit 
observables. By setting the coarse-lattice spacing to be small enough, we ensure that the 
observables are rich enough to extract all the QCD physics. 

Using this blocking framework, I make the following two-pronged argument. I first derive 
the non-controversial claim that exact taste symmetry is recovered in the continuum limit 
of the ordinary, local staggered theory. The concrete physical properties on which this 
conclusion rests are itemized. Like many fundamental properties, no rigorous proofs of 
these physical properties are available. Yet, they are more than just plausible. It is difficult 
to imagine how they could be spoiled without grossly affecting the continuum limit of lattice 
QCD as we know it. 

I next examine to what extent the same physical properties remain valid in the fourth-root 
theory. Properties that have to do purely with the short-distance behavior, and generalize 
relatively easily to the fourth-root theory, include power-counting renormalizability [ill , [l^ . 



131], and the locality of the contributions to the gauge- field effective action generated by the 
integration over the ultra-violet fermion modes. 

A crucial ingredient is the scaling of taste-breaking effects. Because of the absence of 
a local fermion action, it is unclear whether a scaling analysis is at all possible. In other 
words, it is unclear what is the merit of the observation that the taste-breaking effects arise, 
formally, from an irrelevant operator. I bypass this difficulty by building a representation of 
the RG-blocked theory where all the taste violations can be traced back to a local operator, 
whose scahng can be computed by developing the appropriate perturbative expansion. I 
argue that the result should indeed reproduce the scaling law of an irrelevant operator. 

If the taste violations vanish in the limit of infinitely many blocking steps, then, after 
sufficiently many blocking steps, it should be possible to find local theories in the correct 
universality class that provide a good approximation of the fourth-root theory. I construct 
such local theories explicitly via a reweighting of the blocked fourth-root theory that amounts 
to discarding all the taste-breaking terms. Starting from the reweighted theories and working 
back towards the fourth-root theory, I then conclude that exact taste symmetry is recovered 
in the continuum limit of the fourth-root staggered theory as well. As I have explained 
above, this implies the validity of the fourth-root recipe. 

This paper is organized as follows. The RG-blocking framework is introduced in Sec. [ITl 
Technical details have been kept to a minimum, and are mostly relegated to several appen- 
dices. The beginning of Sec. [TTl contains a brief summary both of the section itself and of the 
content of the appendices. Sec. [TTl ends with the introduction of the reweighted theories, at 
which point I give a description of the key steps of the argument, the details of which will 
be presented in subsequent sections. 

The recovery of exact taste symmetry in the continuum limit of the ordinary, local stag- 
gered theory is discussed in Sec. IIIII The reweighted theories are discussed in detail in 
Sec. [IVl Finally, in Sec. El the reweighted theories are used to establish the recovery of 
exact taste symmetry in the continuum limit of the fourth-root theory. My conclusions are 
given in Sec. IVII 

This paper is long, and addresses both conceptual and technical questions. In order to 
help the reader find his/her way through, I have organized the paper such that the essentials 
are summarized in the following (sub) sections: Sec. Ill El explains the technical layout of the 
argument; Sec. IIVDI gives a summary of key properties of the reweighted theories; and 



3 



Sec. I VI I contains the conclusions. For a brief account of this work, see Ref. 
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II. THE RG BLOCKING FRAMEWORK 

In this section I introduce the RG blocking framework, referring mostly to the ordinary 
staggered theory for pedagogical reasons. This section provides a bird's eye overview. It 
contains only the minimum needed to follow the logic of the main argument, as given in 
subsequent sections. Much of what goes into the construction has been relegated to several 
appendices. In the appropriate places, I refer to the relevant appendices for a more elaborate 
discussion. The summary below serves as a "Table of Contents" both for this section and 
for the appendices. 

The blocking transformations are introduced in Sec. Ill Al which also serves to set the 
notation. The first blocking step is special, since it is used to make the transition from 
the standard one-component formalism of staggered fermions to a taste-basis representa- 
tion. This special step is discussed in App. lA II The same appendix also provides technical 
details on the fermion blocking kernels for subsequent blocking steps (App. IA2p . and con- 
tains a proof of the positivity of fermion determinants encountered in the blocking process 
(App. IA3|) . Next, Sec. Ill Bl casts the partition function of the resulting blocked theory in 
a form that will be repeatedly used below. Sec. Ill CI introduces the pull-back mapping of 
operators from the coarse lattice back to the original fine lattice, as well as its main uses. 
A more elaborated discussion of the pull-back mapping may be found in App. [Bl alongside 
with some details on the gauge-field blocking kernels. Also relegated to appendices are the 
generation of blocked gauge- field ensembles (App. [U]), a general discussion of lattice sym- 
metries under the blocking transformation (App. [D]), as well as a more specific discussion of 
the hypercubic and chiral symmetries (App. [E]). 

The reasons why the interacting fourth-root theory cannot be local for any finite lattice 
spacing have been spelled out in Ref. which also contains further details on various Dirac 
operators encountered within the blocking framework. For completeness, a brief review of 
these reasons is given in App. [Fl The remaining two appendices (App. [G]and App. [H]) deal 
with scaling issues. 

In this section, I continue in Sec. Ill Dl with a summary of the main lessons from the free 
theory j20|. Finally, I introduce the fourth- root theory and the various reweighted theories 
in Sec. Ill E| and give an overview of the argument, to be presented in detail in subsequent 
sections, the conclusion of which is that the fourth-root recipe is valid in the continuum 
limit. 



A. The blocking transformations 

Originally, the partition function of the ordinary staggered theory is 

Z = j VUVxDx exp[-5,(W) - xDstaMx] , (2.1) 

where Dstag(p() = Dstag{,x,y;lA) is the staggered Dirac operator, and x{^)-,x{y) is the stag- 
gered field. The fine-lattice coordinates are denoted x,y, and the lattice spacing is a/. The 
link variables are t/^,^, and the gauge field as a whole will be denoted U = {U^^^}. The 
gauge-field action is Sg{U). Summations over all lattice sites will be suppressed. 
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As already mentioned, I assume that all the quark masses are positive.^ In order to avoid 
unnecessarily cluttered notation I will consider the ordinary and fourth-root theories with a 
single flavor of staggered fermions. The generalization to more than one flavor is obvious. 

I will perform n + 1 blocking steps, labeled as k = 0,1, ... ,n. The first, k = step 
is special; it transforms the staggered field from its usual one-component basis to a taste 
basis, which is then retained in all subsequent blocking steps. The k = step maintains 
the number of fermionic degrees of freedom. It is described in detail in App. \M In the 
subsequent blocking steps, thinning out of all degrees of freedom (fermions and gauge field) 
occurs. I have chosen to block (and thin out) the gauge field in the special k = step as 
well, essentially for no better reason than making the notation more tractable. 

In every blocking step the lattice spacing is increased by a factor of two. Thus, = 
2^~^^af for k = 0,...,n. When speaking of the coarse-lattice theory I will refer to the 
theory obtained at the last, k = n step. Its lattice spacing will also be denoted Oc = a„. 
When I increase the number of blocking steps, the coarse-lattice spacing Oc will be held 
fixed in physical units, and the fine-lattice spacing will decrease as = 2~^~^ac. The bare 
parameters on the fine lattice are adjusted to maintain constant physics. I will assume that 
the (fixed) length of each dimension of the coarse lattice is finite. Since, by assumption, all 
quarks are massive, no subtlety should arise in taking the thermodynamical limit. 

For k = 0, . . . ,n, the blocked fermion and anti-fermion fields on the k^^ lattice will be 

denoted as ip^i{x^''^) and ip^ii^^^^) respectively, where x^^^ are the coordinates on the k^^ 
lattice. The indices a and i, both ranging from one to four, are the Dirac and the taste 
index respectively. The blocked link variables will be denoted V^^^'*(fe)- The fc^'^-step blocked 

gauge field ClS db whole is denoted V^^^ = {V^%^}. 

Each blocking step is performed by multiplying the integrand of the partition function of 
the previous step by one, written in a sophisticated form 2^1 (for reviews of the renormal- 



ization group, see Refs. |23l. |24||). Since integrations are over a compact space or else involve 
Grassmann variables, the order of integrations can be chosen at will. The result of the first, 
special blocking step, and the subsequent n ordinary blocking steps, is summarized by the 
following equation: 

Z = J VUVxVx eM-SM-xDstaMx] (2.2a) 

/n 

fc=0 

W^^^r-^^")!)^^"^ exp [- S„(v("\V^("),V^^"^)] , (2.2b) 



where Sn is the final coarse-lattice action. In Eq. fl2.2ap . K, represents the sum of all the 
blocking kernels. The notation {V^'^-', . . .} signifies dependence on the listed fields for all 
< A; < n. Itemizing the blocking kernels. 



A;=0 



ICf + ) , (2.3) 



^ See Refs. [itI . [i9| for a discussion of how to implement the physical theory with a negative quark-mass 



using the fourth-root recipe. 
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in which the subscripts B, F refer to bosons {i.e. the gauge field) and fermions respectively, 
we have 



where 



exp 



and 



(0) 



exp 



/c 



(fc) 



(2.4a) 
(2.4b) 



(2.5a) 
(2.5b) 

(2.6a) 
(2.6b) 



In Eqs. (]2.4bp . (12.5bp and (12.6bp . the fc-range is 1 < k < n. The fermion blocking kernel 
depends on the gauge field on the {k — lattice only. It is an ultra-local and 
gauge covariant linear transformation that maps each 2^ hypercube of the [k — 1)**^ lattice to 
a single site of the k^^ lattice. The blocking parameter ak has mass-dimension one. In this 
paper I will assume that ak is chosen to be of order Note that the gaussian 

integral 

(2.7) 



P^('=W^exp(-4^)) 



yields a trivial constant; in Eq. (I2.2ap . this constant was absorbed into the definition of the 
Grassmann measure. For more details on the fermion blocking kernels, see App. El For the 
gauge-boson blocking kernels, see App. [Bl 



B. A convenient representation 

Keeping track of the taste-symmetry violations is evidently important. For this purpose, 
the action 5"^ that results from the blocking transformations is not very useful because it 
contains multi-fermion interactions. In the fourth-root theory, that action would furthermore 
be non-local, and completely intractable. The problem can be circumvented by noting that 
the fermion blocking transformations are gaussian, as can be seen from Eqs. f l2.2al) and 
(12. 6p . Let us integrate out the original fermion variables as well as all the blocked fermion 
variables, except for those that live on the coarsest lattice. But let us not integrate out the 
original gauge field nor any of the blocked gauge fields. All the integrations we explicitly do 
are, thus, gaussian. The result is 

n „ 

VU\{ \dV^^^] S„(i;W,{V(^)}) / P^("W"^xp (- , (2.8) 
fc=o 
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where 



B, 



n 



(2.9) 



k=Q 



X det' 



k=l 

In the fourth-root theory we will have 



In the ordinary staggered theory one has 7 
|. Here 
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-1 



a. 



0^ + g(°)/^;lg 



(0)t 

(fc)t 



(2.10a) 
(2.10b) 



D 



stag 



(o)tn(o) 



,n. 



(2.11a) 
(2.11b) 



Equation 02.91) includes the constant resulting from the Grassmann integration fl2.7p . for 
each k. The different powers of au arise because the {k + 1)*^ lattice has sixteen times fewer 
fermionic degrees of freedom compared to the fc*^ lattice, except for the = step which 
does not reduce the number of fermionic degrees of freedom. 

For fixed values of the original as well as all the blocked gauge fields, we have the factor- 
ization formula ^22l] 



det(L> 



stag ) 



det(Dn) det 



JJdet 



fc=i 



a 



1/16 



(2.12) 



Equation (12.121) shows how the fermionic short-distance fluctuations are gradually removed 
from the theory. Each factor of det((j^^) results from integrating out fermionic degrees of 
freedom during the fc-th blocking step, and generates an effective action for the collection of 
gauge fields W,V(°),--- M^~'^\ 



cO 
ok 



log det (ao Go) , 
log detfay^*^ Gk 



k = h 



n . 



(2.13a) 
(2.13b) 



All the long-distance physics is contained in the RG-blocked Dirac operator D„. That 
Dn faithfully reproduces the long-distance physics can be seen as follows. By successive 
applications of Eq. fl2.10p . the blocked-field propagator can be expressed in terms of 
the original propagator D'^^lg between a special type of smeared sources built up from the 
product of the fermion blocking kernels (cf. Eq. (1A6P ). As I explain below, this is but a 
special case of a more general mechanism. 

For m > 0, one can show that det{Dk) "-^-'--^i" „„j j„^-/r^-i^ 

App. I A 31 for the proof. 



is strictly positive and det(G';, ) is positive. See 



C. Pull-back mapping 

The RG-blocking transformations start off at the cutoff scale and proceed gradually 
towards lower energy scales. But the blocking transformations facilitate another operation 
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that works in the reverse direction. Suppose that we want to calculate the expectation value 

of an operator O^""^ = O^^'^ {V''"'\ ) defined explicitly in terms of the fields of the n*^ 

lattice. Using Eq. (12. 2p . the expectation value may be calculated as follows. We begin by 
integrating over the n*'^-lattice fields, then over the (n — l)*'^-lattice fields, and so on. If, 
however, we stop at any intermediate step k, the result of the integrations we have done so 
far will be expressible solely in terms of the fields of the k^^ lattice. This procedure defines 
a pull-back mapping of any operator from the n*^ to the k^^ lattice, that by construction 
preserves expectation values. 

Explicitly, for any — 1 < j < n — 1, the pull-back mapping T^^'"-^ : O^"'^ O^^^ is defined 

by 



/n 



exp 



k=j+i 



k=j+i 



C»("^ (2.14) 



As promised, by construction the pull-back mapping preserves the value of observables, 

(^q-i3,n)(^in)^^ ^ (C^"^>„ ■ (2.15) 

Here the expectation value (■■■)„ is defined by the representation of the partition function in 
Eq. (I2.2ap . Taken together, these equations merely say that we may perform the integrations 
in Eq. (12.21) by first integrating over the blocked fields labeled hj j + 1 < k < n, and then 
integrating over the remaining blocked fields as well as over the original fields. The value 
j = —1 accounts for the original fine- lattice theory, and T*^~^'"^ is the pull-back from the 
last-step coarse lattice all the way to the original staggered theory on the fine lattice. 

The pull-back mapping is ultra-local if and only if the blocking kernels are. An operator 
supported on a compact subset of the n}^ lattice is mapped by T^^'"-^ to an operator supported 
on a corresponding, only somewhat bigger, subset of the j^^ lattice. 

An immediate corollary is that the coarse-lattice observables form a proper subset of the 
fine-lattice observables. The coarse-lattice expectation value of the operator O^"'^ is equal 
to the fine- lattice expectation value of the operator T^~^''^''0^'^\ As alluded to earlier, the 
reconstruction of the blocked fermion propagator D^^ from its predecessors is in fact an 
example of the pull-back mapping in action. 

This innocuous corollary leads to another, all important, result. Every coarse- lattice 
observable, being simultaneously a fine-lattice observable via the pull-back mapping, is con- 
strained by all the fine-lattice symmetries. In this sense, the (physical) consequences of the 
exact lattice symmetries cannot "be lost" by the blocking process. 

More can be said on the role of specific lattice symmetries within the blocking framework. 
The interested reader is referred to App. [B] for a more detailed discussion of the pull-back 
mapping. Effects of the blocking transformations on the fine-lattice symmetries, including 
the relevance of the pull-back mapping, are discussed in App. [D] and App. [El 



D. Lessons from the free theory 

In this subsection I review the main results of RG-blocking in the free theory [2^ . With 
no gauge fields, Eq. (12.121) takes a somewhat special form. By a unitary change of variables 
one can switch back and forth between the Dirac operators Dgtag, in the one-component 
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basis, and D taste-, in the taste basis jl, [6|, and^ 



dei{D,tag) = det(Aaste) = det(D„) J] det Gu) '] • (2.16) 

k=l 

An explicit expression for the free RG-blocked Dirac operator Dn may be written down. 
Its taste-violating part A„ (see Eq. f l2.2ip below) has a norm bounded by^ 

||aeA„|| =0(2-") =0(aoK). (2.17) 

In the limit n — > oo, all the taste- violating terms go to zero, and 

lim = Drg^l, (2.18a) 



n— >oo 



lim ^^^(Dtaste) ^ aet\Drg). (2.18b) 

™nLidet(a-^/^^G,^) 

Again, 1 is the identity matrix in taste space. The "one-taste" operator D^g is local, and 
det{Drg) qualifies as a fourth root of det{Dtaste) = det{Dstag) in the sense of Eq. ( I2.18bl) . 



By repeatedly integrating out the short-distance fiuctuations we thus obtain a coarse- 
lattice operator with an exact four-fold degeneracy in the limit n oo. The power-law 
scaling of the taste-breaking terms is clearly as dictated by their origin: irrelevant operators 
with mass-dimension equal to five. Intuitively this can be understood as follows. The 
dimension-five taste- violating terms in Dtaste are multiplied by an explicit factor of oq, the 
initial taste-basis lattice spacing. But the momentum fiowing through the fermion line is in 
effect of order |p| ~ at most. The relative size of the taste- violating terms is therefore 
at most of order ao/oc. 

As a simple corollary of the rigorous work of Ref. 2^1, one can prove that Dk, 



and its inverse Gfc, are all local, bounded operators. Mathematical rigor set aside, one can 
understand how develops an 0{ak) gap directly from Eq. (12. lip . Since the fermion 
mass is very small in any lattice units, it can be ignored for this purpose. Also, from now 
on, I assume that uq in Eq. fl2.10ap has a finite, 0(aQ ^) value. Hence Eq. (12.120 (and not 
Eq. fl2.16p ) must be used also in the free theory. For further details on the /c = step see 
Ref. [15] and App. \M 

For > 1, the massless Dirac operator -Dfc_i satisfies a Ginsparg- Wilson (GW) relation 
26| . in which [75 ^5] takes the role usually played by 75. Here ^5 is the representation of 



75 that acts on the taste index (see App. The eigenvalues of Dk_i thus lie on a circle in 
the right half of the complex plane, with the imaginary axis tangent to the eigenvalue circle 
on the left [l^, lO]. In order to obtain G~^^ from Dk-i, we add the blocking-kernel term 



Q{k)\Q{k) ^ This new term is positive semi-definite. It affects mostly the small-momentum 



Equation (|2.12p reduces to Eq. (|2.16p in the limit ao 00, where det((Q;oGo)" 



1 and D[ 







Dtas 



15|. What I denote as Gk in this paper was denoted Ffe in Ref. 20|, see Eq. (6) therein. The notation 



Gk has a related, but different meaning in Ref. [20]. Also, for compatibility with Ref. [25|, in Ref. 20| the 
blocking parameters ak were all chosen to have a fixed, 0(1) value in units of the coarse-lattice spacing, 
whereas here I make the more natural assimptkm ak — 0(aj7^). 
^ This bound is rigorous in the free theory [20|, [25 1 . 

In the massless limit D^g (Eq. (|2.18p ) satisfies the usual Ginsparg- Wilson relation. 
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modes located near the origin of the Brillouin zone, and pushes their eigenvalues to the 
right by an amount proportional to au- The result is that no eigenvalue remains in an 
0{ak) neighborhood of the origin. In other words, G^^ has developed a gap of order a^-^ 
Its inverse Gk will thus have a decay rate of order au as well. Furthermore, since one may 
also obtain the blocked Dirac operator as 

Z}fc = afc-«^g('=)Gfcg('=)t^ (2.19) 

it follows that the decay rates of the kernels of and of G^]^^ should be 0{ak) too. The 
argument may now be repeated for the [k + 1)**^ step. 

E. Overview of the argument 

When no roots are taken, the local lattice theory with Nf staggered fields belongs to the 
universality class of QCD with 4Nf quarks. The usual universality classification is, however, 
inapplicable to the fourth-root theory, because of its non-locality. There is evidence from 
staggered chiral perturbation theory that the outcome of taking the fourth root may be 
described in terms of an extended Hilbert space containing unphysical states with non-zero 
taste charges: the contributions of these taste-charged states as intermediate states are thus 
also unphysical. A unitary subspace with the correct number of quarks, one per staggered 
field, will exist only in the continuum limit, and only provided that exact taste symmetry 
is recovered [H, [3, M, H 



The range of the non-locality present in the fourth-root theory appears to be set by 
infra-red scales of the theory: the masses of the various staggered pions (isl . 2l| . The taste- 



breaking terms driving the non-locality are lattice artifacts, and would naively be expected to 
vanish in the continuum limit. If indeed exact taste symmetry is recovered in the continuum 
limit of the fourth-root theory, then, when the lattice spacing has become small enough, 
it is logically necessary that local lattice theories in the desired universality class exist 
which provide a good approximation of the fourth-root theory. Such local theories could 
be constructed by simply discarding the taste-violating terms from the blocked fourth-root 
theory. This is the idea behind the introduction of reweighted theories. 

Consider first the ordinary, local staggered theory. With the help of Eqs. (12. 8p and (12.91) . 
the original staggered partition function in Eq. (12. ip can be re-expressed as 

/n 
VU J] [PV(')] Bn (l; W, {V^'^)}) det{Dn) . (2.20) 
fc=0 

In short, this expression results from n + 1 blocking steps, after which the fermion fields 
have been integrated out altogether, while retaining explicitly the integral over the original, 
fine-lattice gauge field W, as well as over the gauge fields V'-'^-* of all the blocking steps. Recall 
that the special = step facilitates the transition from the usual one-component staggered 
basis to a taste basis. 

In order to keep track of taste-symmetry violations, let us split the blocked Dirac operator 



When the eigenvalues Xi are complex we may define the gap as min | A; 
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Dn into its taste-invariant and non-invariant parts: 

Dri -Dinv^n ~1~ 7 (2.21a) 
^inv.n ^inv^n ® 1 5 (2.21b) 

Dinv,n = ^ ^T^ts{Dn) , (2.21c) 

where trts denotes tracing over the taste index only. The taste non-invariant part A„ is 
traceless on the taste index. 

By construction, accounts for physics over distances on the order of the coarse-lattice 
spacing Oc or longer. In particular, its taste violating part A„ accounts for all taste-symmetry 
violations at the energy scale and below. The question is how big are the taste-symmetry 
violating effects in the spectrum of 

Originally, the staggered Dirac operator Dgtag exhibits taste-symmetry violations on all 
scales. In the low-lying spectrum they are small [29]; but they grow gradually with the 
energy scale, until they become 0(1) at the scale of the original lattice cutoff 
I will argue that, nevertheless, by choosing n large enough the taste- violating effects in 
the spectrum of can be made arbitrarily small. RG blocking removes the ultra-violet 
fermionic modes. Their remnant is, mutatis mutandis, the effective action S^jj- (Eq. (12.131) ). 
This effective action is a sum of (products of) Wilson loops and, I will claim, it is local 
on both the ordinary and fourth-root staggered ensembles. Being a local functional of the 
(original and blocked) gauge fields, but not a functional of the fermion fields, it cannot give 
rise to any taste-symmetry violations at large distances. In this sense, the ultra-violet taste 
violations have been eliminated. I will further claim that, with every additional blocking 
step, the remaining taste-violations in the entire eigenvalue spectrum of the blocked Dirac 
operator get smaller, uniformly, basically because this spectrum consists of only "low-energy" 
modes with respect to the fine-lattice scale. In the limit of infinitely many blocking steps, 
taste symmetry is fully recovered. 

We may discard all the taste-breaking effects from the staggered theory, by hand, after 
only a finite number of blocking steps. Truncating the blocked Dirac operator in Eq. fl2.20p 
to its taste-invariant part -DmD,n gives rise to the following reweighted theory 

/n 
VU\{ [W^'^)] S„(l;W,{V('=)})det(An.,n). (2.22a) 

A more conventional looking path-integral representation may be obtained by rewriting 

det(Z)j„j, „) as a path integral over (four-taste) coarse-lattice fermion fields as in 

Eq. (12. Sp . and then integrating out the "tower" of gauge fields except for the coarse-lattice 
gauge field V*-"'^ This gives 

Z,^,^^ = j exp [- 5,„,,„(V("),V'("\^^"^)] , (2.22b) 

which is to be compared with the path integral representation (I2.2bl) of the blocked staggered 
theory. Unlike the staggered theory, the reweighted theory has no shift invariance [3, §]. 
Instead, it has exact taste-f/(4) invariance by construction. Another difference is that the 
above constructed reweighted theory does not have an exact chiral symmetry in the massless 
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limit. ^ 

Reweighting at blocking level n generates a sequence of theories Zinv,n which are different 
from each other, as well as from the staggered theory. But I will argue in Sec. IIIII that, 
because A„ is an irrelevant operator, the (sequence of) reweighted theories has the same 
continuum limit as the (blocked) staggered theory. Each reweighted theory enjoys exact 
taste symmetry by construction, and this implies (the uncontroversial result) that exact 
taste symmetry is recovered in the continuum limit of the ordinary staggered theory. The 
proof works by establishing the existence of a convergent expansion relating the staggered 
and reweighted theories when n is large enough. One must require that all the quark masses 
be non-zero, consistent with the fact that the chiral and continuum limits do not always 
commute [3, [13, H El . 



Moving on to the fourth-root theory, its partition function cannot be represented as an 
ordinary path integral with a local fermion action. Rather, it is given by [10(] 

^root ^ j -pUexpi-S,) det'/^D^tag) , (2.23) 

where the positive fourth root is taken. As in the ordinary staggered theory, this may be 
re-expressed in an n-step RG-blocked form as 

/n 
VU JJ [dV^'^)] Bn (\-M, {V^'^)}) det'/'(Z^„) . (2.24) 
fe=0 

Again let us remove the taste-breaking terms by hand, which gives rise to a new family of 
reweighted theories 

/n 
VU\[ \dV^'^] B^{\-U. {VW}) det(5_,n) . (2.25) 
fc=0 

Here I have used the exact taste symmetry of Dinv,n = Dmv.n ® 1 to take the analytic fourth 
root of its determinant. ^ 

One can represent det(Z)j„„ ^i) as a fermion path integral. This suggests that the validity 
of the continuum limit of the fourth-root theory could be established by, once again, showing 
that the sequence of reweighted theories has the same continuum limit as the blocked fourth- 
root theory. But we must now face two hurdles that were not encountered in the local, 
ordinary staggered theory. 

The new hurdles are addressed in Sec. IIVI First, we must show that the reweighted 
theories derived from the fourth-root theory are local. This is done in Sec. lIVBl by showing 
that, on the basis of plausible assumptions, the effective action S^j^j- and the blocked Dirac 
operator Dk are local on the k^^ lattice scale, on both the ordinary and the fourth-root 
staggered ensembles. Because Dinv,n is defined by a trace projection, Dinv,n and A„ are then 
separately local. 

Introducing coarse-lattice Dirac fields that, this time, carry no taste index, and 

once again integrating out the "tower" of gauge fields except for the coarse-lattice gauge field 



^ More sophisticated reweighted theories may be constructed. See Ref. 15| for a construction that maintains 
the exact chiral symmetry of the m — > hmit. 
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For large enough n, det(-Di„^^„) and dei{Dinv,n) are positive, see Sec. IVl 
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V(") the reweighted fourth- root partition function then takes the form (compare Eq. fl2.22bl) ). 

/n 

VUll [PV(^)] S„(i;W,{V«}) (2.26a) 
k=Q 

^/DV<».pWexp[-S-(v'»',,'"'.5'"')]. (2.26b) 

The "one-taste" action S'™°*„ = 5'[°°*„(V*^"'-*, g^"-*, g'-"'') is comphcated, and contains many 
multi-fermion interactions, just hke Sn and Sinv,n encountered earher. What matters, how- 
ever, is that S'[°°*„ too is local on the coarse-lattice scale if, in particular, Dinv,n is; the 
"surgery" of removing the taste violations has also removed the non-localities of the blocked 
fourth-root theory! 

I stress that the argument for locality of S^jj- and D/^ does not require that the underlying 
theory be local. This non-perturbative argument is very general, and only makes mild use 
of the renormalizability of all theories including the fourth-root theory [ll|, 13] to establish 
the existence of a weak-coupling regime. 

At this point we expect that the reweighted theories Z^°°^^ derived from the fourth-root 
theory are local, and belong to the desired universality class. Now comes the second hurdle. 
Convergence of the reweighted and the staggered theories to the same continuum limit 
depends on the scaling of the taste-breaking effects. But the fourth-root theory does not 
have a local action in the first place, so how are we to perform any scaling analysis? 

The scaling of taste-breaking effects in gauge-invariant observables can, of course, be 
studied numerically, and what one finds is in agreement with what one would naively expect 



30| . But, while no undesirable effects have been encountered at presently accessible 
values of the lattice spacing and the quark masses, numerical results alone cannot alleviate 
the concern that closer to the continuum and/or the chiral limit the scaling of taste-breaking 
effects might eventually be altered in undesirable ways due to the lack of a local lattice action. 

On the lattice, a scaling analysis rests on two pillars. On the non-perturbative side, we 
define a theory, construct local operators within it, and set up an RG transformation. On 
the perturbative side, if the theory is power-counting renormalizable, we can compute the 
scaling of any local operator. An important special case is to take the local operator to be 
the action itself, or individual terms within it. 

In this paper I show how to generalize the scaling analysis to the fourth-root theory. 
First, the fourth- root theory is renormalizable. The multi-gauge- field representation of the 
blocked fourth- root theory, Eq. fl2.24p . then allows us to bypass the lack of a local fermion 
action. Instead, we may study the scaling of the blocked Dirac operator Dk and the effective 
action Having first shown by non-perturbative considerations that both of them are 

local operators, their scaling can be computed by setting up the appropriate perturbative 
expansion, which is done in Sec. IIV CI I find that any local operator scales in the same 
way in the staggered theory and in the corresponding reweighted theories (with or without 
the fourth root). In particular, the taste-breaking part of the blocked Dirac operator, A„, 
indeed scales as an irrelevant operator. 

Finally, in Sec. |V] I reconstruct the rooted theory from the corresponding reweighted 
theories, and establish the validity of its continuum limit. 
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My conclusion rests on renormalizability of the fourth-root theory, concerning which there 
is httle doubt, and the fact that renormahzabihty is "inherited" by the reweighted theories 
(Sec. IIV A|) . My conclusion also rests on two additional key features that have to do with 
locality (Sec. IIV Bp . and scaling (Sec. lIVCp . I give plausible arguments for each of them, 
but confirmation must await more detailed future investigations. Where, then, do we stand 
today? One can draw an indirect but important lesson from the ordinary staggered theory. 
In that case we (believe we) know what is the continuum limit. Moreover, the reweighted 
theories are tightly constrained by the convergent expansion relating them to the - local - 
staggered theory. This leaves little doubt that all the key properties are valid in this case. 
But the argumentation of Sec. II VI makes very little discrimination between the ordinary and 
the fourth-root cases. As I explain in more detail later, this increases confidence that nothing 
essential has been overlooked, and that the claimed properties are valid in the fourth-root 
case as well. 



III. CONTINUUM LIMIT OF THE ORDINARY STAGGERED THEORY 

In this section I discuss the continuum limit of the ordinary staggered- fermion theory. 
Continuum-limit observables that can be computed within the coarse-lattice theory are 
introduced in Sec. IIII Al In Sec. IIIIB| I list scaling relations that follow from a standard 
RG analysis in the ordinary staggered theory. These scaling relations imply the recovery 
of taste symmetry in the continuum limit. This is inferred in Sec. IIII CI by comparing the 
blocked staggered theory to the reweighted theory at each blocking level n. Provided that 
the renormalized quark mass is non-zero, I show that the two theories are connected by a 
convergent expansion when n is large enough, and that any difference between them vanishes 
in the limit n oo. 



A. Continuum-limit observables 

The continuum limit corresponds to the limit of infinitely many blocking steps, which 
is taken while holding fixed the coarse-lattice spacing ac in physical units. The fine-lattice 
spacing a/ of the staggered theory goes to zero, a/A = 2~"~^acA — 0, where A is the 
QCD scale. Constant physics is maintained by adjusting the bare parameters such that the 
renormalized gauge coupling gride) and quark mass mr{ac) are kept fixed. The coarse-lattice 
spacing plays the role of the renormalization scale. The fermion's wave-function renormal- 
ization is controlled by the parameter z^''^ in Eq. (lASp . that fixes the overall normalization 
of the gauge- covariant blocking kernel Q*^^-* in the interacting theory. I will assume that the 
z^'^^^s have been adjusted so as to impose a wave-function renormalization condition on the 
blocked fermion fields at the renormalization scale ac- 

For simplicity I will restrict the discussion to meson observables. Sources for mesons are 
added into the fermion action on the n^^ lattice as follows: 

Si±.M) = i^^^^J■S^-^^P^''\ (3.1) 



For general sources, see appendix B of Ref. 11 1 
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where 

^(„) ^ ^ ^ J^^^in)^ of , (3.2) 



The kernels iSj-"^ are gauge- covariant and uhra- local. Augmenting the fermion action in 
Eq. (12. 8p by the source term (13. 2p and performing the Grassmann integration we obtain the 
partition function with sources (compare Eq. (I2.20p ): 

/n 
VU\{ [PV('=)] B^(l-U, {V('')}) det(D„ + J ■ , (3.4) 

A:=0 

as well as the normalized version 

Zn{ac\ J) = Zniac] J)/Zniac; 0) . (3.5) 

Meson correlation functions, renormalized at the scale Qc, are generated by functional dif- 
ferentiation of ZnicLc', J). The (assumed) existence of the continuum limit means that the 
n ^ oo limit of the normalized generating functional Zn{ac] J) is smooth: 

Zooi^c] J) = lim Z„{ac; J) . (3.6) 

n— >oo 

By differentiation of Zoo{ac', J) one generates continuum-limit meson correlators in euclidean 
space. 

Before moving on let me comment on the set of coarse-lattice observables. In Eq. (12.151) . 
which states the equality of observables under the pull-back mapping, let us choose j = n — 1. 
Instead of using the number of blocking steps as a label, we may use the corresponding lattice 
spacing for this purpose. The equation then takes the form 

(r(a,/2, a,) 0(a,))^^/2 = ((^(«c))„, , (3.7) 

in self-explanatory notation. Equation (13. 7p remains valid in the continuum limit, where 
the pull-back mapping T{ac/2,ac) remains well defined. This equation identifies the ob- 
servables of the coarse-lattice theory Zoo (etc; J) with a proper subset of the observables of 
Zoo{ac/2; J), the coarse-lattice theory with half the lattice spacing. Because observables are 
lost in the blocking, we need to take OcA to be small enough in the first place to ensure 
that the observables derived from Zao{(ic] J) are rich enough to extract all the QCD physics. 
Additional reasons for choosing OcA small will be encountered in Sec. IIV[ 



B. Scaling of irrelevant operators 

I now list the scaling laws needed to establish the recovery of exact taste symmetry in 
the continuum limit of the ordinary staggered theory, as they apply within the RG-blocking 
framework of this paper. 



^"^ See App. lA 21 for the renormalization of composite coarse-lattice operators. 
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First, the very existence of the continuum hmit in QCD derives from asymptotic free- 
dom/^ or, in other words, from the scahng properties of the running couphng QricLc) as a 
function of the bare couphng. 

The remaining scahng laws pertain to the fermion sector. For the restoration of taste 
symmetry in all observables, two scaling laws will be necessary: 

11^" 'II ^ 

A„ < 4 = • 3.9 

These scaling laws are assumed to hold in an ensemble average sense; they do not hold on 
all configurations, but they (are assumed to) hold after averaging over the configurations 
in the ensemble. Each configuration is generated as described in App. O one complete 
configuration consists of a "mother" configuration Ui of the fine-lattice gauge field, as well 
as of "daughter" configurations vl^\ vf^'*, . . . , vf""* of blocked gauge fields. 

The bound fl3.8p relates the lowest eigenvalues of the blocked Dirac operator D„ to the 
running quark mass mr{ac)- This inequality says that the value of mr{ac) measured on an 
ensemble of configurations is set by the lowest eigenvalues of D^, provided that the wave- 
function renormalization condition was imposed at the coarse-lattice scale. The bound (13. 8p 
is needed to tame the infra-red behavior in the fermion sector. We will see in the next 
subsection how it enters. For some further comments on the bound (13. 8p and its parallel in 
the fourth-root theory, see App. [Gl 

The crucial scaling law (13. 9p determines the scaling of the taste-breaking part of 
the blocked Dirac operator Dn (see Eq. (12.211) ). In the taste representation, the leading 
taste violations arise from dimension-five irrelevant operators [1, [g, 20], and the right-hand 



side of Eq. (13. 9p gives the anticipated scaling based on this engineering dimension (compare 
Eq. (I2.17P for the free theory). The "<" sign means that in the interacting theory the 
inequality holds up to logarithmic corrections, that is, powers of log(ac/a/) or, equivalently, 
powers of n. 

In a standard RG application the fermions and the gauge field are both blocked at each 
step, and the scaling laws apply to (the parameters of) the blocked action. Here, the scaling 
laws are assumed within the representation (12. Sp . which is superficially different in that 
none of the gauge fields have been integrated over explicitly. However, because the sources 
couple to coarse-lattice fields only (cf. Eq. (13.30 ). we may imagine that we first integrate over 
W, V^^\ . . . , V*-""^^ and only later over the coarse-lattice (fermions and) gauge field. This is 
equivalent to inserting the source term into Eq. (]2.2bp which contains the blocked-lattice 
action S'„. Thus, in the context of the ordinary staggered theory, the above scaling laws are 
on the same footing as the scaling laws used in a completely standard RG-blocking context. 

Because of the scaling of the taste-breaking part as given by Eq. (13. 9p . the bound (13. 8p 
practically applies to both D:^^ and Dj^^^. When using Eq. (13.80 I will disregard the 
difference between and \\Dl 



mv,n I 



The argument may be extended to other theories such as lattice QED with staggered fermions, if a finite 
but "beyond the Planck scale" lattice cutoff is acceptable. 
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C. Recovery of taste symmetry in the continuum limit 



Assuming the existence of the continuum hmit and the scahng laws (13.81) and (13.91) . I 
will now prove that exact taste symmetry is recovered in this limit for all the coarse-lattice 
observables of the ordinary staggered theory, provided mr(ac) > 0. The proof makes use of 
the reweighted theories introduced in Sec. IIIE[ and reveals why it is necessary to avoid an 
exactly massless lattice theory. 

I first add the source term (13. ip to the n*'^ reweighted theory (cf. Eq. (I2.22ap ): 

/n 
VU\[ [dV^'^)] S„(l;W,{V(^)})det(An.,n + (3.10) 

I also introduce a family of partition functions^^ (with sources) in which D„ is replaced by 
Dinv,n + tAn, whcrc t takcs values in the interval < t < 1. Explicitly, 



/I L 
VU\{ [dV^^^] B„(i;W,{V('=)}) (3.11a) 

X det [d^ + J- 5^")) det 1 + (t - 1) A„ {d^ + J ■ S^""^^ 

/n 
I^W JJ [W^'^)] B„(^1;W,{V('=)}) (3.11b) 

fc=0 

X det (An.,n + J ■ det 1 + tAn (An.,n + J " ^^^^j 

These partition functions interpolate between the reweighted theory for t = 0, and the 
staggered theory for t = 1. Normalized generating functionals are defined in analogy with 
Eq. (13. 5p . and the n oo limits in analogy with Eq. (13. 6p . 

I will now show that the n — > cxd limit does not depend on t, viz., 



^ooiflci J) Winter, cxi(t y J) ^inv,oo(^(^cj J) • 

The key step is to bound the last factor in Eq. (13.1 lap as 

exptr log (1 + (t - l)A„D„i) = 1 + ©(e^) , 

where 



D 



-1 



A, 



(3.12) 



(3.13) 



(3.14) 



I have used that A„ is traceless on the taste index, as well as the geometric series for the 
staggered propagator = D~^^^ — D~^^ „ A„_D~^ n + ' ' ' • The sources are infinitesimal and 
do not interfere with any bound valid for J = 0. The (unnormalized) partition function of 
any interpolating theory, Z inter, n(t,ac; J), is obtained from the staggered partition function 
by reweighting with the left-hand side of Eq. (13.130 . Because the scaling laws (13. 8p and (13. 9p 
hold by assumption on the staggered ensemble, it follows that 



2-(n+l) 



(3.15) 



For positivity of the determinants see Sec. |Vl 
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We arrive at several important conclusions. First, each term in the Taylor expansion of the 
logarithm in Eq. fl3.13p is bounded by the corresponding power of e„. Second, because the 
product actrir^ac) is held fixed, there will be an tlq such that, for any n > uq, one has < 1 
and the Taylor expansion converges. It follows that the change in any (meson) observable 
over the interval < t < 1 is 0(aj). Finally, since e„ — > for n — > oo, we readily arrive at 
Eq. dsn. 



The reweighted theories Z,, 



inv,n V"^C) 



J) have exact f/(4) taste symmetry by construction. 



and the same is true for the limiting theory ^mD,oo(a.c; J)- But the limit is independent of t, 
and so the staggered generating functional Zoo{clc, J) has exact taste symmetry as well; the 
RG-blocked staggered theory becomes taste-invariant in the limit of infinitely many blocking 
steps if, in particular, Eq. (13.91) holds. As far as the rate of restoration of taste symmetry is 
concerned, this is recognized as the familiar result that the discretization errors of staggered 
fermions are proportional to the (fine) lattice-spacing squared jl, 31, IH, 33, 34 . 

In the massless staggered theory one has mr(ac) = 0, and the bound (13.151) becomes an 
empty statement. Therefore it is not possible to infer the recovery of full taste symmetry in 
the exactly-massless case. This is consistent with the established fact that the continuum 
and the chiral limits of staggered fermions do not always commute 16|, llTl, llSl, [19 . 

The result I have established readily generalizes to all the coarse-lattice observables. 

Given a coarse- lattice operator = C(")(V'(''\ V^")), I introduce the notation {of^^ 

for a "mixed," unnormalized expectation value where the sea quarks have Dirac operator 
Dint,n{t) = Dinv,n + whilc the valence quarks have Dirac operator Dint^nit'). To be pre- 
cise, the Boltzmann weight is given by Eq. (13. lip (with J = 0), while the valence propagators 
are given by D~^^^{t'). We then have 



(n) 



Cj"^ exp 



Of^)^ {l + 0{el)) 



trlog (l + (t - l)A„Z}-i 



(3.16a) 
(3.16b) 



Ot^) (l + 0(e„)) 



(3.16c) 



Equation (I3.16al) re-expresses an unnormalized observable of Z inter, n{t) as a correlation func- 
tion computed on the staggered ensemble with a reweighting factor. Equation (]3.16bp follows 
from the validity of the scaling laws on the staggered ensemble. Last, Eq. (I3.16cp replaces any 
valence propagator Dj^^ „(t') by the staggered propagator D~^. The expansion of D^^^ ,^{t') 
as a power series in AnD~^ has the same convergence properties as that of the logarithm in 
Eq. (I3.13p .^^ As a special case, (1)^ = (1)^^ (1 -fO(e^)), and Eq. (13.160 readily generalizes to 
normalized expectation values. It follows that all the staggered and the reweighted (t = 0) 
observables have the same n — > oo limit, which now establishes the exact taste symmetry of 
all observables. 

Equality of all the observables implies the equality of the "fixed-point" coarse-lattice ac- 
tions obtained in the limit n — oo. Comparing once again the staggered and the reweighted 



I have allowed for an 0(e„) mismatch in the observable, though presumably for any physical quantity of 
interest one could construct a coarse-lattice observable that would have only 0{aj) discretization errors, 
in which case the additional mismatch incurred in Eq. (|3.16cp is likely to stay O(e^) as well. 
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theories, this means (cf. Eqs. fl2:2b|) and fl2.22bp ) 



S^{ac;V,ip,i!) = \im Sn{ac;V,ip,i!) 

= Urn Sinv,n{ac]V,ip,ij) ^ Sinv,cx>{ac;V,ij,ij) . (3.17) 

n— >oo 

In analogy with Eq. (13.71) . I have dropped the blocking-step label attached to the coarse- 
lattice fields, and traded it with the coarse-lattice spacing. Adding in the source term fl3.ip 
then gives 

Zoo{ac;J) = J VWi^Vi, exp [-Soo (ac; V, ^,^) - ■ 5(a,; V) V^] . (3.18) 

By Eq. f l3.17p . the action S'oo(ac; V, ■?/', ■?/') of the limiting RG-blocked staggered theory has 
acquired exact taste-?7(4) invariance. Equation (13.180 can be used to derive taste-5't/(4) 
Ward-Takahashi identities that will be exactly satisfied in the limiting theory. The corre- 
sponding result in the fourth-root theor y ( see discussion below Eq. (15.151) ) will put on a 
firm basis the observation made in Refs. |lll . [l9j that no paradoxes can be derived from the 



extended taste symmetry of the continuum-limit fourth-root theory. 



IV. REWEIGHTED THEORIES AT WEAK COUPLING 

In the previous section I have shown that, for large enough blocking level n, the reweighted 
theory can be reached from the local staggered theory by means of a convergent expansion. 
The converse is also true: Setting t = 1 in Eq. (]3.11bl) . we can reconstruct the staggered 
theory from the reweighted theory. The convergence of the (i:-)expansion in Eq. (13.11bp is 
controlled by ej„^,^„, where 



inv,n 



(4.1) 



The notation || • ||j„^ means that the norms are now to be evaluated on the reweighted 
ensemble. In view of the established scaling of A„ on the staggered ensemble, Eq. (13.91) . 
the difference between any staggered-ensemble expectation value and the corresponding 
reweighted-ensemble expectation value must be very small. Indeed we must have einv,n ~ ^n, 
up to corrections which are of higher order (in either of them). 

The ability to go back and forth between the staggered and reweighted theories implies 
that the reweighted theory associated with the local staggered theory must have the following 
key properties: 

• A suitable notion of renormalizability; 

• Locality on the coarse-lattice scale; 

• Validity of the scaling laws of Sec. IIII Bl including in particular the scaling of A„ as 
in Eq. (13. 9p . on the reweighted ensemble. 

In this section I explore direct evidence for these key physical properties. The main output 
is that, step by step, every argument about the (four taste) reweighted theory derived 
from the local staggered theory generalizes straightforwardly to the (one taste) reweighted 
version of the fourth-root theory. In Sec. |Vl this will allow me to establish the validity of 
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the continuum limit of the fourth-root theory. The n ^ oo hmit of the blocked fourth-root 
theory will be reached via the corresponding limit of the sequence of reweighted theories. 
Since the reweighted theories are all local on the coarse- lattice scale, the same will be true 
for the (common) limiting theory. 

Power-counting renormalizability of the ordinary and fourth-root staggered theories, 
alongside with the derived reweighted theories, is discussed in Sec. IIV A[ Locality of the 
reweighted theories is addressed in Sec. IIVBI The scaling laws are discussed in Sec. IIV C\ 
relegating some further issues to App. [Hi I summarize the emerging physical picture in 
Sec. UVDl 

In the rest of this paper I will assume that the coarse-lattice scale has been chosen to 
satisfy -C A~^. This has the following implications. (1) Because of asymptotic freedom, 
the running coupling constant gr{a-c) is weak at the coarse-lattice scale as well as on all 
shorter distance scales. (2) One can define lattice-regularized QCD to be local if it is local 
at the coarse-lattice scale. (3) The coarse-lattice observables are rich enough to extract all 
of the QCD physics. 



A. Renormalizability 

I begin with a brief account of what is known about the renormalizability of the ordi- 
nary and fourth-root staggered theories (Sec. IIV A ip . I then offer a natural definition of 
renormalizability for reweighted theories, from which it follows that a reweighted theory is 
automatically renormalizable if the underlying staggered theory is (Sec. IIV A2|) . 



1. Staggered theory 



Unlike Wilson fermions [35j, the task of deriving lattice power- counting theorems and all- 



orders renormalizability remains to be completed for staggered fermions (for recent progress. 



see Refs. [13|, l36|). Still, it is widely believed that the ordinary staggered theory is renor- 
malizable to all orders. The main evidence comes from a one-loop calculation accompanied 
by the observation that the staggered-fermion symmetries forbid the generation of any rel- 
evant or marginal terms not already present in the staggered action [8[]. In particular, the 
taste-breaking terms remain irrelevant to all orders. 



As first noted in Ref. 12 



all-orders renormalizability should extend from the ordinary 
staggered theory to the fourth-root theory. The argument relies on the familiar replica 
trick. Consider an ordinary staggered theory with copies, or replicas, of equal- mass 
staggered fields. At this stage, Ur is a positive integer. At each order in perturbation theory, 
the counter-terms needed to renormalize the lattice theory will be polynomials in n^, because 
Ur only enters as an overall multiplicative factor attached to every closed fermion loop. Next, 
we consider the analytic continuation in Ur to arbitrary real values, which corresponds to 
raising the staggered determinant to a (possibly fractional) power n^. This continuation is 
unique, because the (polynomial!) n^-dependence of the diagrammatic expansion at each 
order, including the counter-terms, is already known. Thus, the counter-terms derived for 
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The relevance of the rephea trick [33| for the low-energy pion sector of the fourth-root theory has been 
discussed in Refs. [ISl, [27 
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integer rir will be just enough to renormalize the fractional-power theory for any value of 
Ur- (While this captures the essence of the argument, it amounts to an over- simplification. 



For a more thorough discussion, see Ref. 



Thus, while the fourth-root theory is non-local [15|, renormalizability is not lost. Re- 
taining renormalizability turns out to be the absolutely essential starting point from which 
everything else follows. 



2. Reweighted theory 

Renormalizability of a lattice theory means that, by adjusting the bare parameters, the 
correlation functions of the renormalized fields have a finite limit when the lattice spacing 
goes to zero and momenta in physical units are kept fixed, to any order in perturbation 
theory. In the present RG-blocking context I will assume that the correlation functions 
under study are constructed from coarse-lattice fields. The external momenta all belong to 
the Brillouin zone of the coarse lattice. 

While the reweighted theories depend on both the fine- and coarse-lattice scales, I will 
adopt exactly the same criterion to define when they are renormalizable. The implication is 
that renormalizability carries over automatically from the staggered theory to the derived, 
reweighted theories. The reason is that, to leading order (in the fine-lattice spacing), the 
difference between a given coarse-lattice correlation function in the staggered and in the 
reweighted theory is equal to the taste-breaking part the staggered correlation function. In 
Sec. IIV CI below I will argue that the taste-breaking part of any diagram vanishes in the 
continuum (n oo) limit, and so the staggered and the reweighted theories assign the same 
continuum-limit value for every correlation function,^'' to any order in perturbation theory. 
This prediction applies to both the ordinary and the fourth-root staggered theories, in fact 
to any real value of the number of replicas n^. 

In order to avoid unrelated non-perturbative complications, as well as to ensure the 
existence of a weak-coupling regime, I will further restrict Ur to positive values where the 
one-loop beta function (which depends linearly on n^) remains asymptotically free. The 
main results of Sec. IIV Bl below are valid under these mild restrictions only. As a preparation 
for the next stage, let me write down the staggered-fermion RG-blocked partition function 
(without sources) for a general number of replicas n^. 



/n 
VU J] [w^'^)] Bn (n.; W, {V^"')}) det"'- (A™,™ + tA„,) , (4.2) 
fc=0 



in which I have also kept the interpolating parameter < t < 1 of Eq. (13.111) . All the 
partition functions studied in this paper are special cases of ZJ^^ Ur, ac). As claimed above, 
they are all renormalizable. 



The notion of renormalizability is often assumed to include the requirement that the continuum-limit 
value of each diagram can be made equal to the value computed using some other regularization method 



by a finite renormalization [35| . To the extent that this is true for the staggered theory, this will readily 
generalize to the (sequence of) reweighted theories as well. 
Extending this claim to the non-perturbative level is the subject of Sec. |Vl 
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B. Locality of SK^ and 



The main result of this subsection is that both the effective action S^j^j (Eq. (12.131) ) 
and the blocked Dirac operator D^. (Eq. (I2.iup ) are local on the k^^ lattice scale. In more 
detail, integrating out fermionic degrees of freedom at the k^^ blocking step generates a local 
effective action S^jj- for the gauge field, and the Dirac operator governing the remaining 

fermionic degrees of freedom is local too. Because -DmD,n is defined by trace projection, 
Dinv,n and A„ are separately local. I will argue that this is true on the ensemble generated 
by Zn(t,nr,ac) defined above, for any t and (in the indicated ranges). The argument 
relies on the renormalizability of the lattice theory defined by Z„(t, nr,ac), but it does not 
require that that theory be local in itself. 

As noted in Sec. Ill El (see Eq. (12.261) ). a corollary of crucial importance is that the theory 
defined by Zn{t, rir, ac) turns out to have a local coarse-lattice action whenever raising of the 
fermion determinant to the n^J^ power is an analytic operation. A local coarse-lattice action 
defines when a reweighted or an interpolated theory is local. A local coarse-lattice action 
is obtained for the ordinary (ra^ = 1) staggered theory, as well as for any theory derived 
from it by varying t; and for the = 1/4 theory at t = 0, which is recognized as the (n^^) 
reweighted theory derived from the fourth- root theory, cf. Eq. (12.261) .^^ 

Let me begin with RG-blocking in a pure Yang-Mills theory. While again no rigorous 
proofs exist, it is widely accepted (see e.g. Ref. [l^]) that the pure Yang-Mills lattice action 
obtained after n blocking steps is local. The coarse-lattice action can be approximated by 
the continuum form 

Sn{0'c) ~ -IT, — T / F^u + discretization errors. (4.3) 
gr{(^c) J 

The ensemble of coarse- lattice configurations generated by the Boltzmann weight Bn{a^ = 
exp{—Sn{ac)) will correspond to the correct running coupling gr{ac)- (See App. Ofor the 
generation of blocked-lattice gauge field configurations.) 

Next let us consider lattice QCD. After integrating over all fields except for the coarse- 
lattice gauge field, one arrives at a Boltzmann weight of the general form (again using 
continuum notation) 

Bn{ac) ^ exp 

Here "0(1)" stands for terms occurring at zeroth or higher order in the expansion in powers 
of Qrio-c), which are also of zeroth order in the expansion in powers of the lattice spacing 
{i.e. terms that survive the continuum limit). The form (14.41) depends on renormalizability, 
because the ultra-violet divergent part of the fermion determinant that renormalizes the 
coupling constant has been separated out explicitly. All other effects of the integration over 



The same is true at t = for rir = Us/'^, where Us is a positive integer, interpreted as the number of 
equal-mass sea quarks. All other (fractional!) values of Ur fail to yield a local action, even at t = 0, 
because there is no local Dirac operator D such that det^''"''(D) = det(I?i„i,^„). Notice that, in Sec. Illli 
it was not necessary to use the locality of the reweighted theories because the ordinary staggered theory 
is by itself local. 



9'^ 



d X F.,,, + 0(1) + discretization errors 



(4.4) 
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the fermion fields are contained in the last two terms in Eq. fl4.4p . These terms obviously 
include the effects of virtual quark loops, and they are not local. For the considerations 
in this subsection, the only thing that matters is that the renormalized Yang-Mills action 
inside of the Boltzmann weight is parametrically larger by 1/ g'^{ai^. Equation fl4.4l) applies 
to the partition function well (after all but the integration over V*^"-' has been 

done), because this partition function too defines a renormalizable theory. 

With Eq. fl4.4l) in hand, we are ready to discuss the locality properties of S^j^ and D^- 
Let me summarize the relevant discussion from Sec.[TTl(in particular around Eq. (12.191) ). but 
now in terms of the hermitian [l^ operator 

Hk = [l,®i,]G-^\ (4.5) 

First, functional differentiation of S^jj with respect to the (original or blocked) link variables 
generates expressions that depend on both and H^^ . Locality of S^j^ then follows 
provided that and if^^ are both local on the k^^ lattice scale. The locality of Dk-, Hk 
and Hj^ is established iteratively using Eqs. (]2.11bp and (I2.19p . (For the A; = step, see 



Ref. [15|].) The only non-trivial step is to demonstrate the short-range nature of Hj^ ^. In 
the free theory, this follows because has an 0{ak) gap. 

In the interacting theory one has to replace the notion of a spectral gap by the notion of 
a mobility edge. The properties established in the free theory will carry over to any smooth 
gauge field, and, more technically, to any order in lattice perturbation theory. This leaves 
open the following question. In the presence of very rough, lattice-size structures in the 
gauge field, or "dislocations," could Hk develop much smaller eigenvalues, which in turn 
would spoil the short-range character of its inverse? I will now argue that the answer is 
negative. 

Before coming to the main argument I should note that it is logically possible that all 
eigenvalues A^-'^'' of may always satisfy {Xf^'^l > A^]„ > 0, for some A^]„ = 0{ak)- In 
other words, might have an 0{ak) gap for all gauge fields. If this is true, we are done. 
Because it is unknown if this is true, I will disregard this possibility. 



According to the theory of disordered systems (see Ref. |39|, |40] and references therein), 
the right question becomes what is the mobility edge of Hk- In general, the spectrum of 

(k)i (k) 

Hk will consist of localized eigenstates (at the scale a^) with eigenvalues < |A. | < Ac . 
Above the critical value Xi'^"' the eigenstates become extended. By definition, Xi^'^ > is the 
mobility edge of Hk- The value of Xi^^ is a property of the ensemble. 

I will argue that the mobility edge of Hk is 0{ak) in the weak-coupling regime. On general 
grounds, any localized eigenmodes lying below the mobility edge, be their eigenvalues as 
small as they may, will not spoil the short-range character of the inverse H-^ |3g,|4l|. The 
decay length of Hj^^ is thus 0{ak) as required. 

At scales where lattice QCD is weakly coupled, the physics that goes into the mobility 
edge is simple. Because it may be unfamiliar, and since Hk itself has not been studied 
numerically yet, let me digress to describe results obtained in the study of the mobihty edge 



of the hermitian Wilson operator Hy\/. In Ref. j40|] the mobility edge of Hw was determined 
for the super-critical bare mass arriQ = —1.5, and for a range of values of /3 = G/qq on 
quenched ensembles. This example is relevant for the following reason. First, for the chosen 
parameters, Hw and Hk both have 0{a~^) gaps in the free theory, where a is the relevant 
lattice spacing. Moreover, at any super-critical bare mass the spectrum of Hw can, and does, 
reach zero. Therefore, the analogy will become relevant in case that future numerical work 
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will demonstrate the existence of low-lying eigenvalues of Hk. According to the argument 
below, the corresponding eigenmodes will necessarily be localized. 

For the case at hand, the most interesting finding of Ref. was that the mobility edge 



was very close to the free-theory gap for several different gauge actions, even at the not- 
very-large cutoff scale ~ 2 GeV. At stronger coupling (lower cutoffs) the mobility edge 
did go down, eventually reaching zero when the Aoki phase was entered. 

Returning to the RG-blocking context, we set the coarse-lattice spacing such that gricLc) 
is as weak as we like. Therefore we are interested in values of the mobility edge at weaker 
couplings than any of those already studied. The results of Ref. jiO] suggest that, on a very 
weakly-coupled pure Yang-Mills ensemble, the mobility edge of any operator, including the 
super-critical Hw and Hk, will be very close to the free-theory gap; for f] ^ oo the mobility 
edge will continuously approach the free-theory gap. 

It remains to consider the inclusion of a fermion determinant raised to some positive 
(but not necessarily integer!) power at weak coupling, as described by Eq. f l4.4l) . This 
should have little effect on mobility edges which are already 0{a~^). A different power of a 
fermion determinant does change the beta-function and the running of the gauge coupling, 
so we should be a bit careful in what we are comparing. Consider first the operator if„ 
of the last blocking step. We may compare the value of its mobility edge on a pure Yang- 
Mills ensemble to the corresponding value on a dynamical-fermion ensemble that has the 
same coarse-lattice running coupling gr{ac)- Once the coarse-lattice couplings are equal, the 
remaining contribution of the fermion determinant in Eq. (14.41) is parametrically smaller by 
a factor of g1{ac) compared to the Yang-Mills action which is the leading term. Changes 
to the spectrum near a mobility edge which is already 0{a^^) and, therefore, any further 
deviations of the mobility edge itself from the free-theory gap, are expected to be very small. 

We actually need to know something about the mobility edges of Hk for all k < n. For 
each k, we may compare the blocked staggered ensemble to a new Yang-Mills ensemble 
chosen such that the running coupling gr{ak) at the k^^ lattice scale is the same in the two 
theories. Again a similar conclusion will follow. 

What can, and will, be significantly affected by the inclusion of fermion determinants 
is the small-eigenvalue localized spectrum (if there were any near-zero eigenvalues to begin 
with; see for example Ref. [s^). Since the Boltzmann weight contains det(G^^) = det{Hk) 
raised to a positive power, the transition from the pure Yang-Mills ensemble to dynamical 
staggered- fermion ensembles will lead to fewer near- zero eigenvalues of Hk, for all k. 

Let me summarize the anticipated physical situation. In the free theory, Hk has an 0{ak) 
gap. On weakly-coupled pure Yang-Mills ensembles, the mobility edge of Hk is expected to 
be very close to the free-theory gap and, thus, 0(0;^) by itself. Now, starting at rir = 0, 
let us gradually increase the power of the staggered- fermion determinant in Eq. (14. 2p . while 
maintaining a fixed renormalized coupling gr{cLc) at the coarse-lattice scale. Any near- zero 
eigenvalues of Hk will be gradually suppressed, but otherwise nothing much should change 
in the spectrum of Hk, for all k. In particular, any further change in the mobility edge will 
be even smaller than the, by itself small, deviation from the free-theory gap on the pure 
Yang-Mills ensemble. As a result, the mobility edge of Hk will remain 0{ak), and the 
decay length of Hj^^ will remain 0{ak)- As explained above, this implies the locality of Dk, 
Hk, H^^, and S^jj, for all k < n 



I use asymptotic freedom to bound gridk) by gr{ac)- This restricts my conclusions to the range of 
rir-values where the one-loop beta-function is negative. See also footnote [T3l 



24 



C. Scaling of A„ 

In Sec. IIV Al I have explained why the partition function (14. 2p defines a renormahzable 
theory for any and t. For generic values of these parameters, this theory is non-local. 
Nevertheless, thanks to its power-counting renormalizability we may study the scaling of any 
local operator constructed within such a theory. This includes in particular the operators 
listed at the end of the previous subsection, whose locality I have just established by non- 
perturbative considerations. 

I will be mostly interested in the scaling of A„, the taste-breaking part of the blocked 
Dirac operator, in the staggered (t = 1) and in the reweighted {t = 0) theory. In this 
subsection I argue that Eq. (13.91) correctly describes the leading power-law scaling of A„. I 
furthermore find that the logarithmic corrections to the scaling of A„ depend on Ur but not 
on t. The argument is heuristic, and will have to be confirmed by future calculations. 

Setting up perturbation theory is in principle straightforward. In reality, perturbation 
theory for a reweighted theory is unfamiliar, and technically rather different from ordinary 
staggered perturbation theory. It can be gradually built up in several steps: 

Step 1. Ordinary staggered perturbation theory; 

Step 2\ Staggered perturbation theory with the fermions in a taste basis ob- 



tained via a unitary change of variables p, |32 



Step 2. Staggered perturbation theory with the fermions in a taste basis ob- 



tained via a gaussian smearing RG-like step [15|; 

Step 3. Multi-gauge-field perturbation theory for the blocked staggered theo- 
ries of Eqs. ([21]), fl2:20|) or ^TIM : 

Step 4- Multi-gauge-field perturbation theory for a reweighted theory (or, more 
generally, for Z„(t, n^, Oc)). 

As I will explain, computing the scaling of A„ within the fully developed perturbative setup 
of Step 4 can be reduced, via Steps 3 and 2 (skipping Step 2'), to a calculation in ordinary 
staggered perturbation theory (Step 1). Focusing on the fourth-root theory, the outcome is 
that in spite of the lack of a local action, we nevertheless have at our disposal a local operator 
A„ that, on the one hand, accounts for all the taste violations in blocked observables, and, 
on the other hand, is controlled by staggered perturbation theory. 

Ordinary staggered perturbation theory is a well-developed technique and there is no 
need to discuss it here (see e.g. Ref. [lH and references therein). I now discuss all the other 
steps listed above. 

In taste-basis perturbation theory the fermion momentum ranges over a reduced Brillouin 
zone, and the sixteen sites of each 2^ hypercube are accounted for by the Dirac and the 
taste degrees of freedom. A taste-basis perturbative expansion is usually not used for two 
reasons. The unitary transformation to a taste basis is not unique, and complicates the form 
of many symmetries jl, [32| (see Apps. El [D] and [E] for more details). Also, in a taste-basis 
diagrammatic expansion, taste violations occur in both the propagator and the vertices, as 
can immediately be seen by inspection of the free taste-basis Dirac operator (see App. [Hlfor 
further discussion). This is to be contrasted with the usual staggered perturbation theory, 
where the momentum-space propagator is taste-symmetric and taste violations reside in the 



vertices only [8|, [11 
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Next, consider the staggered theory obtained by performing the special k = "blocking" 
step introduced in Sec. Ill Al for the fermions only. No blocking is applied to the gauge field U. 
The resulting Dirac operator is Dq of Eq. fl2.10al) . Following the notation of Ref. pj| where 
its explicit form was derived, in this subsection I denote it as D^aste = -Dtaste(«0) "^)- In the 
free theory one can write Dfaste = i{A + B) + Ai, where A, B and M. are all hermitian, and 
correspond to the Dirac ® taste structures 7^^®!, 75®^^^^5, and 1®1 respectively. As usual, 
DtasteioiQifn) satisfies a GW relation in the massless limit. A key feature of Dtaste{otQ)^) 
is that, provided ao < oo, dropping the taste-breaking part B does not introduce any new 
doublers into the theory [l5| . The gaussian-smearing transformation reduces to the previous 
unitary transformation in the limit ao ~* oo, where M.{oo,m) = m. Here, as in Sec. Ull I 
assume that ao = 0{a~j^). 

With its extra technical hassles, staggered perturbation theory with the Dirac operator 
Dtaste is, clearly, highly relevant to the present blocking framework. For the purpose of 
the discussion below, I only need to draw attention to one fact. By using the general 
procedure of Ref. [i^ one can prove that Dtaste retains all the staggered symmetries, albeit 
in a complicated form.^^ These symmetries forbid the appearance of taste-breaking relevant 
or marginal terms through loop corrections jl] . The taste- violating part of the fermion self- 
energy, denoted Ft.^., is therefore O(p^) in lattice units. But since taste violations occur now 
in both vertices and propagators, this result necessarily represent a delicate cancellation. 
Schematically: V^^ — ^ = OijP') where both V^^ and F'^ are 0(1), and where FJ'^ 
accounts for the contribution of diagrams with at least one taste-breaking vertex, while F' ^ 
accounts for the remaining contributions, in which taste-breaking arises from the fermion 
lines only. 

The thing to notice is that this delicate cancellation would be hampered had we truncated 
the free propagator {-i{A + B) + M) / {A^ + B"^ + M"^) to the "linearized" form {-i{A + B) + 
M.)/{A? + M."^) obtained keeping only the first taste-breaking term in the geometric-series 
expansion of the free propagator about the taste-symmetric {—lA + M.)/{A^ + M."^). This 
truncation is not entirely foolish, because it does not introduce doublers into the theory. Yet, 
it mutilates shift symmetry (see App.[H]for a related discussion). Dynamically, the reason is 
that, for lattice-scale momenta that contribute to the loop integrals, the taste-breaking part 
of the propagator is not small relative to the taste symmetric part. Therefore, truncating 
the propagator will give rise to large, imbalanced changes in V^^ and F'^ . 

I will now argue that the main change brought about by (many!) iterations of the RG 
transformation is that the taste-breaking term A„ becomes uniformly small, for all the 
coarse-lattice momenta. As a result, for large n the scaling of taste-breaking effects will be 
controlled by diagrams with a single insertion of A„. 

Before we do any RG-blocked diagrammatic calculation, we must first set up the appro- 
priate perturbative expansion. Usually lattice perturbation theory is based on the expansion 
of the link variables as U^^^ = exp(i(7aA^^a;). When the definition of the lattice theory makes 
use of the representation fl2.8l) . a similar expansion will have to be applied to the entire 
tower of gauge fields W, V'-^-', V*-^\ . . . , V*-""-* simultaneously. With this, one can in principle 



The underlying reason is that, on top of the previously applied unitary transformation, gaussian smearing 
changes the propagator by a contact term only. The long-distance propagator is unchanged, and therefore 
all the symmetry constraints on long-distance correlation functions must hold, as can be shown using 
the pull-back mapping (Sec. Ill Cl and App. iPjl . The modified, Ginsparg-Wilson-Liischer (GWL) chiral 



symmetry 



43 1 associated with the GW relation is a special case of the general construction of Ref. [42 



26 



set up the perturbative expansion for every theory that can be cast in the form of Eq. (14.21) . 
because the closed- form expressions for Dk, Dk, and S^jj- as functional of all the gauge 
fields are known. 

The next stage is to consider the multi-gauge-field perturbative expansion of the RG- 
blocked staggered theory (t = 1 in Eq. fl4.2l) ). This perturbation theory will reproduce all 
the scaling laws derived using ordinary staggered perturbation theory. The reason is simply 
that, as already noted in Sec. IIIIBl one can first integrate out U, V^'^\V^^\ V^'""^). At 
this point, one has effectively recovered the diagrammatic expansion derived from the action 
Sn (Eq. (]2.2bp ) of the n^^ level blocked theory. The scaling of the parameters of Sn, in 
turn, must agree with the predictions of ordinary lattice perturbation theory. While this 
description is, strictly speaking, applicable for integer values of n^, all other values can be 
reached via the replica trick, cf. Sec. IIV A[ ^^ 

The last step is to show that A„ must scale in the same way in the (blocked) staggered 
theory and in the corresponding reweighted theory. This result is established by isolating 
the diagrams that determine the scaling of A„ in each theory, and showing that they amount 
to exactly the same set of diagrams. 

In the blocked staggered theory, A„ gives rise to taste breaking effects either through 
vertices or though the expansion of the free blocked propagator as = „ — 

D~^^ ^/S.nD^nv n "I — ' • Ncxt comcs the main observation. The momentum that flows though 
any fermion line is, in all cases, a coarse lattice momentum p < l/a^. After sufficiently 
many blocking steps, any coarse-lattice momentum will be very small in fine-lattice units: 
paf ~ af/ttc <C 1. (This is true whether the fermion line forms a closed loop or connects 
to an external leg.) In contrast, A„ embodies taste-breaking effects coming from all higher 
momentum scales up to the fine-lattice cutoff l/ct/. This means that any mechanism needed 
to ensure the smallness of all taste breaking effects on the coarse lattice, such as cancellations 
based on symmetries, must be built into the functional form of A„ = A„(W, V*-"-*, . . . , V*-'^-') 
itself. Said differently, most of the needed cancellations must occur over distance scales 
much smaller than Qc- Therefore, they will not occur in an expectation value with multiple 
insertions of the operator A„, unless they already occurred in every expectation value with 
a single insertion of A„. 

Within the multi-gauge-field diagrammatic expansion, this translates into the statement 
that any insertion of A„ in any diagram must scale as aj /a^ in coarse-lattice units. Given 
a diagram of the blocked staggered theory, let us now drop every contribution where the 
total number of insertions of A„ (coming from both propagators and vertices) is bigger than 
one. In the remaining taste-violating diagrams, the fermion propagator is D~^^^, and they 
contain exactly one insertion of A„. These diagrams determine the scaling of A„ in the 
blocked staggered theory, and, as argued above, will reproduce the taste-breaking scaling 
laws of the original staggered theory. 

But, clearly, the very same set of diagrams is what determines the scaling of A„ in the 
reweighted theory! I conclude that, in any reweighted theory, A„ must scale in the same 



Alternatively, we may consider the pull-back of A„ to the original fine-lattice staggered theory, where we 
may again study its scaling as a function of n. 

As long as we stay within the confines of perturbation theory, this procedure gives meaning to the blocked 
action Sn for any real value of rir- Conceptually, this is similar to the way the diagrammatic expansion 
gives meaning to the dimensionally-continued action in dimensional regularization. 

The argument can be generalized to t =/= 0, with the same conclusion. 
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way as in the original staggered theory. While this will not be needed, the diagrammatic 
correspondence is clearly tight enough to encompass the logarithmic corrections as well. 
Thus, the logarithmic corrections depend on n^, but not on t. 

More generally, the scaling behavior of any local operator must be the same in the stag- 
gered and in the reweighted theory, simply because the scaling will be insensitive to dropping 
A„ from the blocked Dirac operator. (The special case considered above amounts to tak- 
ing the local operator to be A„ itself.) This implies that the physical observables of the 
reweighted theory should have scaling violations proportional to aj, as predicted by (pertur- 
bation theory for) the original staggered theory. The n oo limit therefore yields a "perfect 
action" theory [2J]. (See, however, App. [D]for a discussion of related technical issues.) 

In summary, we learn two important lessons. After many blocking steps, A„ will be small 
on any staggered or reweighted ensemble. We may thus compute its scaling behavior on ei- 
ther ensemble by appealing to perturbation theory for the multi-gauge-field representation 
of the blocked staggered theory {t = 1). Also, as long as we allow for coarse-lattice observ- 
ables only, this calculation further reduces to a conventional scaling calculation in staggered 
perturbation theory, augmented by the replica trick (for non-integer rir). 

In particular, I find that the power-law scaling (13.91) is valid in the reweighted theory 
derived from the fourth-root theory. Notice that I have assumed that A„ scales like a 
dimension-five (and not like a dimension-six) operator. But, as explained in Sec. IIII Ct 
thanks to taste-tracelessness of A„ this assumption is consistent with 0(aj) scaling of the 
taste-violating effects in all the physical observables. In Sec. El the scaling of A„ will be 
used to establish the validity of the continuum limit of the fourth-root theory. 

In this subsection I gave only a very minimal discussion of the multi-gauge-field dia- 
grammatic expansion. In App. [Hi I illustrate some further aspects of this expansion by 
considering a few examples of terms which are expected to occur in A„. 

D. Summary and future work 

We are almost done. In the next section, the reweighted theories will be used to establish 
the validity of the fourth-root theory in the continuum limit. This conclusion is a straight- 
forward corollary of the emerging physical picture of the reweighted theories. I therefore 
pause to summarize what has been learned. 

All the reweighted theories introduced in Sec. Ill El share the following key features: (1) 
renormalizability, (2) locality, and (3) the same scaling laws as the underlying staggered 
theory. At the starting point is the all-essential observation that the fourth-root theory is 
renormalizable if the ordinary staggered theory is. From this point on, basically the same 
reasoning was applied in both the ordinary and fourth-root cases. In fact, for the most part 
the arguments generalize to any real number of replicas within the range specified above 
Eq. g3). 

As briefly discussed in Sec. IIV A It renormalizability of both the ordinary and fourth-root 
staggered theories is not as solidly established as in other cases (notably Wilson fermions). 
But there is no real reason to doubt it either. For a recent, more thorough discussion, see 
Ref. As explained in Sec. IIV A 21 the reweighted theories "inherit" their renormaliz- 

ability from the underlying staggered theory, in a rather trivial way. 

Locality of the reweighted theories at the coarse-lattice scale rests on the locality of the 
effective action 5'^^^ and the blocked Dirac operator Dj., on the relevant ensemble. Those 
locality properties, in turn, are set by the range of H^^, where the hermitian operator Hk 
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(Eq. (14.51) ) accounts for the short- distance fermion modes integrated out at the k^^ step. I 
have argued in Sec. lIVBI that what is needed is that the mobility edge of Hk be 0(1) in units 
of the k^^ lattice scale. I have drawn an analogy to a recent application of the theory of 
localization to lattice QCD, specifically, to a study of the mobility edge of the super-critical 



Wilson operator [39|,|40|. I concluded that, thanks to the existence of a weak-coupling regime 
(which in turn is a consequence of renormalizability) , both the mobility edge of Hj^ and the 
range of Hj^^ will be 0(1) in k^^ lattice units, as required. Obviously, it will be necessary 
to confirm the claims by numerical investigations of Hk itself. The first non-trivial instance 
is provided by the k = 1 blocking step.^^ 

Our knowledge about the ordinary staggered theory strengthens the claims I have made. 
In the ordinary staggered theory, based on standard RG considerations one assumes that the 
blocked action Sn (Eq. fl2.2bp ) will be local on the coarse-lattice scale. For this to be true, 
the locality properties of and must be as claimed. But my reasoning in Sec. IIVBI 
did not discriminate between the ordinary and the fourth-root ensembles. This lends higher 
credibility to the proposed physical picture in the fourth-root well. 

My claims are on stronger footing for rir = as well: this quenched limit is closer to the 
actual setup of the work reported in Ref. The fourth- root value Ur = 1/4 may thus be 



reached by interpolation, starting either from = or = 1- Once again, this supports 
the claims made in the fourth-root case. 

I now turn to the scaling of the taste-breaking effects represented by A„. The basic 
difficulty is simply that, in the fourth-root theory, there is no local fermion action. Thus, it 
is unclear if the taste violations that reside in the fermion sector are amenable to a scaling 
analysis. 

First, the fourth-root theory is renormalizable. Therefore, even though the theory is 
non-local, we have a power counting and we can study the scaling of any local operator. 
Specifically, I have shown that a scaling analysis in the fermion sector is made possible 
thanks to the multi-gauge-field representation introduced in Sec. [Tll This gives us access to 
the operator A„ that accounts for all the taste-symmetry violations in blocked observables. 
According to the discussion of Sec. lIVBt A„ is a local operator; therefore its scaling can be 
computed using the appropriate (multi-gauge-field) perturbative expansion. Finally, I have 
argued that the needed scaling calculation ultimately reduces to a calculation in ordinary 
staggered perturbation theory augmented by the replica trick, and that A„ indeed scales as 
an irrelevant operator should. 

My arguments in Sec. IIVCI were heuristic, and it is clearly necessary to confirm them 
by performing the appropriate perturbative calculations. The actual scaling of A„ can also 
be investigated numerically, at least on the fourth-root ensembles provided by MILC ji^ ]. 
A first study was performed last year Most of the arguments of this section rely on 
being in a (sufficiently) weak-coupling regime, and it is important to understand how close 
are we to this region in practice. Numerically reweighting is clearly a challenge, which, if 
successfully tackled, could further strengthen confidence in the entire framework. Another 
challenging project is to perform an accurate comparison of the predictions of the various 
perturbative expansions to numerical results obtained e.g. by measuring Wilson loops [sgI] 
or by adapting the Schrodinger-functional technique [i^, 46, ll^, 48]. 



Because of special features of the fc = step, the operator Hq is guaranteed to have a gap in the interacting 
theory to t^. T he same is not true for k > I. 
26 See Ref. 



11| for similar examples taken from Condensed Matter physics. 
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V. CONTINUUM LIMIT OF THE FOURTH-ROOT THEORY 



Assuming the properties of the reweighted theories discussed in Sec. llVt in this section 
I prove the vahdity of the fourth- root theory in the continuum hmit. As explained earher, 
when the blocking level n is high enough, one can either reach the reweighted theory from 
the staggered theory via a convergent expansion, or work the other way around. I find it 
appealing to reconstruct the fourth-root theory from the reweighted theory, because the 
latter is local, and is already expected to be in the correct universality class. The argument, 
that otherwise follows the same logic as in Sec. IIIIl is given in Sec. IV A[ I add several 
comments on the scaling analysis in Sec. IV B[ 

A. Recovery of locality in the continuum limit 

In the fourth-root theory I will make use of the scaling laws (compare Eqs. fl3.8p and 

m) 

l|An||_ < 5' (5-1) 

where the subscript "mf " refers to the reweighted ensemble, and 

"''inv,n 

where 

minv,n = rririac) + 0{af/al) . (5.3) 

The scaling of A„ was discussed in Sec. IIVCI The leading power-law behavior of A„ is 
robust. It is unchanged by taking the fourth root, and it is also independent of reweighting. 
Turning to Eq. fl5.3p . the origin of the rightmost term is simply that the transition from 
Dn to Dinv,n amouuts to dropping A„. The latter is 0{af/a1) which, therefore, could entail 
similar changes in the eigenvalues. See App. [G] for some further comments on the bound 
(15. 3p . Similar considerations show that the effective coupling constant of the reweighted 
theory, denoted ginv,n, satisfies 

9L,n = 9riac) [1 + ((a//ae)')] , (5.4) 

where I have used Eq. (14.41) and the taste-tracelessness of A„. This implies that the 
reweighted theory is in a weak-coupling regime if the coarse-lattice staggered theory is, 
and vice versa. 

I will restrict the present discussion to meson observables of the fourth-root theory. The 
generalization to all other observables requires additional technical steps which are discussed 
in Appendix B of Ref. pJJ]. The physical meson observables are taste singlets. They are 
probed by restricting the source term of Sec. IIII Al to the form J ■ iS*^"-*, where now iS^-"^ = 
[^("^ ® 1]. Here ^J"^ carries no taste index and, as usual, 1 is the identity matrix in taste 
space. Switching notation from J to J is meant to remind us that the sources now couple 
to taste singlets only. The blocked fourth-root partition function with these sources is given 
by 

/n 
VU n Sn(i; W, {V«}) det^/^(D„ + J ■ . (5.5) 
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Equation (15. 5p means that observables are constructed as follows [27]. Fermion - anti- 
fermion contractions are done in the same way as in the ordinary staggered theory; then one 
applies the extra "replica" rule that a factor of ^ is to be attached to every closed fermion 
loop occurring in the observable itself (in other words, to every valence staggered-fermion 
loop). With this replica rule in place, the pull-back mapping defined in Eqs. (12.141) and 
(IBip remains valid, and the same is true for Eq. (12.151) . Thus, the ultra- local nature of the 
pull-back mapping is preserved, even though the lattice action itself is not local. 
I add the same source term to the reweighted theories, which gives rise to 

/n 
VU\{ \dV^'^] Bn[l;U, {V^'^}) det(A„.,. + J ■ S^'^^) . (5.6) 
k=0 

Here I have used the exact taste invariance of the reweighted theories and the taste-singlet 
nature of the sources to take the analytic fourth root. In analogy with Sec. Illlt I also 
introduce interpolating theories (with the same source), whose partition functions can be 
expressed as: 



/n 

fc=0 



(5.7) 



xdet (An.,„ + J ■ det^/" 1 + tA„ (A™,™ + J ■ 5^")) 



Normalized varieties of all partition functions are defined in analogy with Eq. (13. 5p . e.g. 

^inv,ni^c'i J) = ^inv,n{'^c] J ) / Z i^v ,ni^c'i 0) . (5-8) 

I will assume that the continuum limit of the (sequence of) reweighted theories exists: 

2Ltoo(«c; J) = lim Z™°*„(a,; J) . (5.9) 

As usual, this is based on the scaling of the coupling constant itself, which, in turn, is only 
negligibly affected by reweighting (cf. Eq. (I5.4p ). 

We are now ready to reconstruct the observables of the fourth-root staggered theory from 
those of the reweighted theory. To this end I use that, on the reweighted ensemble. 



exp 



^tr log(l+tA„Ain) 



l + 0{ite,n.,ny) . (5.10) 



The definition of einv,n is the same as in Eq. (14. ip . except that this is now in the context of 
the fourth-root theory of course. The similarity between Eqs. (I3.13P and (I5.10p is clear. It 
follows that ^ ^ 

Zi:?Ljt,a,;J) = ZXnia,;J) [l + O {{te^Y)] . (5.11) 

Because mr{ac) scales logarithmically while A„ is suppressed by a power of the fine-lattice 
cutoff (cf. Eq. (15. ip ). it is guaranteed that, for n above a certain value, we will have einv,n < 1 
and, with it, convergence of the t-expansion in Eq. (I5.10p . Once again, einv,n — ^ for n — * oo, 
and the continuum limit is independent of t, 
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We now recall that the reweighted theories are local on the coarse-lattice scale, as can 
be seen from the path integral representation (12.261) . and that they belong to the correct 
universahty class. In the hmit n — oo we thus have 

^;ntoo(«c; J) = J VVVqVq exp [sZU^c, V, g, q)-qJ- S{a,; V) q] , (5.13) 

where the limiting action 

SZU^c;V,q,q) = hm SZ'^^ia,;V, q,q) , (5.14) 

is local too. But, by Eq. (15.121) . Zl°f^{ac\ J) accounts for the continuum- limit observables 
of the (blocked) fourth-root theory as well. This establishes the validity of the continuum 
limit of the fourth- root theory. 

1 conclude with two additional observations. My first comment concerns the physical 
consequences of the continuum-limit taste symmetry of the fourth-root theory. One can 
lift the restriction on the sources and consider (meson) observables with a general taste 
structure. The reweighted partition function then takes the form 

/n 
PWn[^^^'^] ^n(i;W,{V('=)})detl/^(An.,n®l + (5.15) 

Taste-5'f/(4) Ward-Takahashi identities may now be derived by varying J. These identities 
will be exact in the reweighted theory for every ra, and will be true up to 0(e^„^ „) corrections 
in the blocked fourth- root theory. In the n ^ oo limit, these identities will become exact 
in the fourth-root theory as well. If, however, we reinstate the restriction to taste-singlet 
observables, then Eq. (I5.15P evidently reduces to Eq. (15. 6p . This means that no paradoxes 
can be derived based on the taste symmetry of the continuum- limit fourth-root theory 
(as claimed in Ref. [49|). The taste non-singlet states live in an extended, non- unitary 
Hilbert space; but a unitary, physical sub-space exists. A more practical conclusion concerns 
deciding when one is allowed to use taste non-singlet operators (such as those in Eq. (lD2p ). 
which are often advantageous numerically, instead of taste-singlets ones. For a detailed 



discussion of these issues, see Refs. 11, 19, 21 



Another observation is that, as 1 have assumed in Sec. Ill El det{Dinv,n) and det(-Dj„t),n) 
will both be strictly positive when e^, einv,n < 1. 1 begin by re- writing the blocked staggered 
determinant as 

dei{Dn) = det(A„.,n) det(l + A^^-^ . (5.16) 



Because [75® ^5] -Dn is hermitian [15||, and Dinv,n = Dinv,n®1- accounts for its taste-invariant 
part, it follows that 75-Dj„^_„ is hermitian too. Therefore, det{Dinv,n) and det{Dinv,n) are 
real. Moreover, by my assumption, the expansion of the determinant in Eq. (15.160 is 
convergent (on both the staggered and reweighted ensembles), and the rightmost deter- 
minant in Eq. (15.160 is thus strictly positive. Since det(D„) is strictly positive too (see 
Sec. IllBp . it follows that det(Dj„t,^„) is strictly positive. Next 1 consider det(i5j„^^„). Be- 
cause det'^{Dinv,n) = det{Dinv,n)^ we know that det{Dinv,n) cannot be zero and cannot flip 



^"^ As explained in Sec. lIVCi S^°°*^ is a "perfect" action. 
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sign. By considering the limit where the bare mass of the original staggered theory goes to 
infinity, it follows that det(Dinv,n) is strictly positive. 

Finally I should note that it is quite certain that the bounds I have made use of in this 
section and in Sec. [Ill] must represent over-estimations. I return to this point in Sec. IVTl 

B. Scaling in the reweighted theories revisited 

A key result of this paper is that the fermion sector of the fourth-root theory becomes 
amenable to a scaling analysis by means of the multi-gauge-field representation of the blocked 
theory. As explained in Sec. IIV C\ the scaling analysis in the fourth-root theory can be carried 
out by reducing it to a calculation in ordinary staggered perturbation theory augmented by 
the replica trick. 

Interestingly, in the reweighted fourth-root theory, the needed scaling laws may be found 
without making any reference to the replica trick in staggered perturbation theory. Ac- 
cording to this alternative route, the calculation of the scaling of A„ (which is still done as 
described in Sec. IIVCP proceeds by first considering only reweighted theories with = Arir 
quark species, where Ug is a multiple of four, and therefore is integer. This means that 
the complete calculation, including the part done on the staggered-theory side, involves 
local theories only. The scaling of A„ for any other number of quark species Ug in the 
reweighted theory can now be found without any further reference to the staggered theory. 
We simply analytically continue the previous result to the desired value of n^. As usual, 
because the dependence is known in closed form, the analytic continuation is uniquely 
determined. Unlike in the staggered theory, however, this analytic continuation only relates 
local (reweighted) theories to other local (reweighted) theories! In particular, the scaling 
in the one-taste reweighted theory is inferred from the scaling in reweighted theories where 
the number of quark species is a multiple of four, without ever having to perform a scaling 
analysis in the fourth-root theory. 

Thus, this line of argument relates the needed scaling properties of the reweighted fourth- 
root theory to the local, ordinary staggered theory, while passing only through local theories 
at intermediate steps. The fourth-root theory is then encountered only at the very last stage, 
where we reconstruct it from the reweighted theory, as was done in Sec. IV A[ 

I comment in passing that, "forgetting" where they came from, the reweighted theories 
Zinv,n or Z^°°^^ each constitute a family of local theories defined on a lattice with spacing 
Oc, which depend on an additional parameter n. The role of this parameter is similar to the 
fifth dimension L5 of domain-wall fermions: when either n or L5 are sent to infinity, a GWL 
chiral symmetry is recovered. The actual construction of the reweighted theories would 
amount to a gross "overkill," if our only aim was to find solutions of the GW relation. The 
merit of the construction is that the same local operator, A„, controls both the violations 
of the GWL chiral symmetry (that originates from the staggered f/(l)e symmetry) in the 
reweighted theory, and the deviations of the latter from the corresponding staggered theory. 



Of course, in a one-flavor theory, a GWL symmetry exists only in the free theory. 
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VI. CONCLUSION 



Like a journey through a dark wood, when deahng with a difficult problem in quantum 
field theory one can never be too sure which is the right way, and where danger is lurking. 
I have concluded that the fourth-root recipe is valid in the continuum limit using plausible 
assumptions. Plausibility is, at the end of the day, in the eye of the beholder. In this 
concluding section, I give my personal perspective on what has been gained. 

In a way, this paper trades one set of questions for another. But while at the starting 
point the questions were rather vague, the new questions are focused, technical, and testable. 
The initial worries basically stem from our lack of experience with non-local theories. A 
formal expansion of the staggered action suggests that the taste-breaking terms are irrelevant 
operators, that would naively be expected to vanish in the continuum limit. But it is unclear 
how to perform a scaling analysis when there is no local fermion action in the first place. 
Related, one must also translate the (tentative) claim "locality is recovered in the continuum 
limit" into a well-defined statement. 

This paper offers a solution to these problems. By first RG-blocking the (fourth-root) 
staggered theory and then enforcing exact taste symmetry by reweighting, we obtain local 
coarse-lattice theories in the desired universality class, which provide a good approximation 
of the (fourth-root) staggered theory once the number of blocking steps is large enough. 
That the reweighted and staggered theories are indeed close to each other, follows from a 
scaling analysis, which, in the fourth-root case, is made possible by the multi-gauge-field 
representation of the blocked theory introduced in Sec. [TTl Within this representation, the 
taste-breaking effects all arise from the taste-breaking part A„ of the local, blocked Dirac 
operator D„. 

The reasoning of this paper has been presented early on in Sec. Ill El and I now recapitulate 
it: All-orders renormalizability of the reweighted theories follows from that of the (ordinary 
and fourth- root) staggered theories (Sec. IIV A|) : making mild use of renormalizability to 
establish the existence of a weak-coupling regime, a robust non-perturbative consideration 
shows that the reweighted theories are local (Sec. IIV Bp : the scaling of the local operator 
A„, which embodies all the taste violations in the blocked theory, can be traced back to a 
calculation in ordinary staggered perturbation theory (augmented by the replica trick in the 
case of the fourth-root theory), and the result is that A„ indeed scales as an irrelevant oper- 
ator (Sec. IIV CI) : the smallness of A„ on the reweighted ensemble enables the reconstruction 
of the staggered theory from the reweighted theory by means of a convergent expansion 
(Sec. |V]); in the continuum limit, the difference between the (blocked) staggered theory and 
the reweighted theory, which is already known to be in the correct universality class, van- 
ishes. For the ordinary, local staggered theory this implies that exact taste symmetry has 
been recovered; for the fourth-root theory, this implies that it has become local. Thus the 
fourth-root theory provides a valid regularization of QCD. 

This conclusion depends on confirming the key properties of the reweighted theories. 
This amounts to verifying their locality, checking the actual predictions of their perturbation 
theory, as well as testing these predictions non-perturbatively (by numerical methods). A 
summary of what each of the above amounts to has been given in Sec. IIVDI 

A detailed-level comprehensive study of all the properties of the reweighted theories would 
be a major endeavor. Nevertheless, already now there is good reason to believe that the 
fourth-root theory is indeed a valid regularization of QCD. This conclusion derives from the 
comparison to the local four-taste staggered theory. In short, our understanding of the local 
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staggered theory is on essentially the same footing as with any other local lattice fermion 
method. If the continuum limit of the four-taste staggered theory is what we think it is, 
then it is difficult to see how the claimed properties of the reweighted theories derived from 
it could go wrong. But then, the arguments for the key properties of the reweighted theories 
apply, basically unchanged, to the four-taste (derived from local staggered) and one-taste 
(derived from fourth-root staggered) cases, both of which constitute local theories. Thus, it 
is also difficult to see how any key property of a reweighted theory could go wrong in the 
one-taste case, if this does not happen in the four-taste case. 

By following this line of argument one can in fact avoid any reference to the replica trick 
- which is the manifestation of the non-locality in staggered perturbation theory. Instead, 
one first derives the scaling laws for reweighted theories where the number of quark species 
is a multiple of four (and, thus, the original staggered theory is local). From this, one infers 
the scaling laws for reweighted theories with any other number of quark species (Sec. IV Bp . 
Thus, any reference to the non-local fourth-root theory is avoided until the very last step 
where it is reconstructed from the reweighted theory. 

Taste-breaking effects in the spectrum of the staggered Dirac operator are largest at 
the (fine-lattice) cutoff scale. But the largest taste-breaking effects are not a major source 
of non-locality; in fact they entail basically no non-locality, because RG blocking trades all 
ultra-violet fermion modes with a local correction to the gauge-field action. This observation 
is nothing but a (part of the standard) description of how symmetries broken by the lattice 
regulator are recovered in the continuum limit. A key result is that this feature is not lost 
by the fourth-root theory. 

The remaining non-local effects have been argued to be associated with the dimensionless, 
small parameter a/ A (isl . 21, 27, 2^. It follows from the results of this paper that all the 
non-local effects should be controlled by (powers of) a/A. In Sec. [Till and Sec. |Vl I have 
bounded the relative size of taste-breaking effects in long-distance observables, hence also 
the relative size of non-local effects, by powers of aj /{almr{ac))- But, because the coarse- 
lattice spacing Oc is basically arbitrary (apart from the restriction ac ^ A~^), this has got 
to be an over-estimation. In all likelihood, the actual relative size of the taste-breaking 
and the non-local effects is on the order of afA'^/niphys (or powers thereof), where rriphys 
is the renormalized quark mass extract from some low-energy observable. This is based 
on the anticipation that, on low-energy modes of the staggered Dirac operator, the actual 
magnitude of taste-splittings among quartets of eigenvalues should scale like a/A^. For 



related theoretical discussions, see Refs. [15|, Il9|, |50l . 

In numerical simulations, the taste-symmetry violations are observed to decrease rapidly 
as a/ is decreased, and indeed to be roughly proportional to (a<i(a/) ci/A)^ (with "improved" 
staggered quarks) j^. The presence of a fixed physical scale, and not some 3> A, makes 
it possible to extrapolate to the continuum limit using present-day computer resources. 

This work lends strong support to the physical picture advocated in Refs. BHiJ: The 
non-localities of the fourth-root theory can be interpreted in terms of an extended Hilbert 
space containing states with, in general, non-zero taste charges. The physical subspace 
consists of the taste-singlet states. The exact taste symmetry, recovered in the continuum 
limit, relates physical and unphysical states, and its Ward-Takahashi identities play a crucial 
role in establishing unitarity in the physical subspace. 

It is interesting to consider a closely related problem, namely the use of the fourth-root 
recipe for finite-density simulations (see Ref. j50| and references therein). Here there is a new. 



three-fold difficulty. First, there are all the general difficulties having to do with a complex 
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measure, that set in when Re (/x) 7^ 0, where fi is the chemical potentiaL Second, when trying 
to apply the fourth-root recipe one confronts phase ambiguities, and the systematic error 
they introduce must be kept under control. Last, a non-zero quark mass is no longer an 
effective infra-red cutoff. The eigenvalues reach the origin in the complex plane for realistic 
values of the chemical potential. Even in this case, it has been argued in Ref. HO] that 



everything is in principle under control, provided that the continuum limit is taken before 
the thermodynamical limit. The crucial grouping of eigenvalues into quartets - near the 
origin in the complex plane and beyond - can still be done when one is close enough to the 
continuum limit, if the volume (in physical units) is finite. However, the systematic error 
due to the phase ambiguities is parametri cally much larger, and grows with a positive power 
of the volume. For more details, see Ref. [50]. 

Returning to zero density, the up and down quark masses used in numerical simulations 
are larger than their physical values. Extrapolation of numerical results to the physical point 
requires the appropriate low-energy effective theory. For the development of staggered chiral 
perturbation theory (S^PT) see Refs. [H, 34, 11]. Recently, based on plausible assumptions 
within the context of the chiral effective theory, it has been argued that S^PT augmented 
by the replica trick is indeed the correct low-energy description of the pion sector of the 
fourth-root theory 

It will be interesting to re-derive S^PT with the replica trick directly from the underlying 
theory, the (RG-blocked) fourth-root theory. The difficulty is that the effective theory 
depends on the number of replicas Ur both explicitly, as well as implicitly through the n^- 
dependence of its low-energy constants. Normally, the dependence of low-energy constants 
on the parameters of the underlying theory is non-perturbative, and is not known. The 
challenge is to cast the (RG-blocked) underlying theory into a new form where the necessary 
analytic continuation in the number of fermion species of some type can be done in a closed 
form. Work on this subject is in progress (sH . 
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This difficulty never appears wlien Re (^) = 0. When the quark mass is strictly positive the staggered 
determinant is strictly positive too, and the positive, analytic fourth root can always be chosen. 
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APPENDIX A: THE FERMION BLOCKING KERNELS 

In this appendix I describe the fermion blocking kernels in some more detail. The tran- 
sition to a taste representation in the special k = step is discussed in App. lA 1[ The 
blocking kernels of all subsequent steps are introduced in App. lA 2[ In App. lA 3[ I prove the 
positivity of det(D„) and det(G~^). 



1. Taste representation in the interacting theory 

In the case of free staggered fermions there is a unitary transformation between the 
one-component field x{^) and the taste-basis field ipai{x^^^), given explicitly by [1,0] 

^Sh^^'^) = E (7[^7?7?7?)«. X(2x(°) + r) . (Al) 

rM=0,l 

For notation see Sec. Ill A[ Writing Eq. (lAip compactly as = F x, the "conjugate" Grass- 
mann variables are related by ^p = xF^ Recall that the fine-lattice spacing a/ of the one- 
component field, and the lattice spacing oq of the taste-basis field, are related via oq = 2a f. 
Equation (lAip makes use of the embedding of the taste-basis lattice into the fine lattice, 
and the fact that each fine- lattice site has a unique representation as = 2x^/^^ + r^, where 

In a free theory, RG blocking normally works by suppressing modes with a lattice-scale 



momentum. It is therefore natural to apply the blocking in the taste basis [20|, where all 
the long-distance physics comes from the vicinity of the origin in the Brillouin zone. The 
one-component formalism would be inconvenient^^ because the long-distance physics comes 
from all sixteen "corners" of the Brillouin zone. 

Unlike the free theory, an equal choice between a one-component basis and a taste basis 
does not exist in the interacting theory. Lattice QCD with staggered fermions must be de- 
fined in the one-component formalism. According to power-counting arguments and explicit 
one-loop calculations, only this formalism has enough symmetry to ensure the multiplicative 
renormalization of the staggered-fermion mass term and the recovery of full rotation and 
taste symmetries in the continuum limit js], [s^] ■ This state of affairs poses a difficulty for the 
RG program. The question is how to accommodate all the symmetries of the standard one- 
component formalism in a taste-basis representation that will, in turn, provide the starting 
point for the succession of RG-blocking steps. 

The unitary transformation from the one-component basis to the taste basis can be 
promoted to a gauge-covariant one, 

i^^h^'"') = (A2) 

^■^=0,1 

Setting Xo = 2x''^^ for short, an explicit choice for the parallel transporter is ^ 
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RG blocking of free one-component staggered fermions was discussed in Ref. 
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where jl is the fine-lattice unit vector in the /i direction. With the notation of Eq. (lA2p we 

similarly have ip^i i^^^'') — xQai\^^^^)^ where hermitian conjugation applies to the color 
matrices. Parallel transporting the fine-lattice variables entails well-defined transformation 
properties for the taste-basis variables under fine-lattice gauge transformations. Thus, gauge 
invariance is maintained by the blocking transformation. 

The covariant blocking kernel flA2p illustrates, however, an inherent problem. Any con- 
crete choice of the gauge- covariant blocking kernel will transform non-trivially under hyper- 
cubic rotations. In Eq. flA2p this is seen both in the special role of the hypercube's site with 
relative coordinates = 0, because these relative coordinates transform non-trivially under 
hypercubic rotations; and, for a similar reason, in the specific ordering of traversing the axes 
in Eq. flA3|) . 

The solution adopted in this paper is to perform the transition from the original one- 
component formalism to a taste representation as a gaussian RG-like transformation, in 
which no thinning out of the fermionic degrees of freedom (but only of the gauge field) 
occurs. For the taste-basis Dirac operator resulting from this transformation, see Ref. (l5| . 
This comprises the special k = blocking step introduced in Sec. Ill A[ Within the gaussian 
blocking transformation, ^(x) is loosely equal to Q^^\x^^^;U) x, and ip{x) is loosely equal 
to X W). For a precise statement, see App. [Bl 

Of course, replacing the unitary change of variables flA2l) by a gaussian transformation 
does not by itself solve the difficulty with hypercubic rotations. But we may now overcome 
it by creating a coherent superposition over a family of different blocking transformations. 
This is explained in App. |E] (see also App. iDl) . 



2. Fermion blocking kernels for 1 < k < n 

For completeness, let me specify the fermion blocking kernels of the subsequent, 1 < k < n 
blocking steps. In the free theory I take 

g(fc)(5.W)^(fe-i) = l y ^{^-i)(2#) + r). (A4) 
16 

r-M=0,l 

In analogy with Eq. ( 1A2I) we may define a covariant version, 

gW(5(fc).V(^-i))V;('^^i) = — J2 w(25('=),2x('=) + r;V('-i))^('=-i)(2x('^) + r), (A5) 

'■m=0,1 

where the parallel transporters are defined analogously to Eq. (1A3|) . but now in terms of 
the blocked gauge field of the {k — 1)*"^ lattice. These definitions imply that the linear 
transformation Q^^^ is unitary, whereas for 1 < k < n, the product Q^^'^Q^^^^ is equal 
to (2;*^'^))^/16 times the identity matrix on the k^^ lattice. The difficulty with hypercubic 
symmetry recurs at every blocking step, and again it is solved in a similar manner (see 
App. El). 

The constants z^^'' are adjusted to impose a wave-function renormalization condition on 
the fermion fields at each blocking level. Usually, lattice renormalization produces factors 
of log(a/i) where /x is the renormalization scale. But in an RG-blocking setup one has 
a flfc-i, /i — l/ofc at the k^^ blocking step. Whence log(a/i) log(2), and we may expect 
= 1 + Cfc5f^(afc)/(167r^), where Ck = 0(1). Of course, the product of all the z^'^^'s can 
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diverge (or vanish) in the hmit n — ^ oo, as dictated by the integrated anomalous dimension 
of the fermion field. 

The expectation value of any product of local composite operators constructed from the 
coarse-lattice fields will always be finite. Therefore, composite operators do not necessarily 
require a separate renormalization. (Once again this can be explained by the fact that 
the ratio of the cutoff and renormalization scales is a finite, fixed number.) One might, 
however, opt to impose specific renormalization conditions for certain composite operators. 
A renormalization condition imposed on a composite operator at the coarse-lattice scale will 
in general entail some finite renormalization. 



3. Positivity of det(L'„) and det{G-^) 

Here I prove that det(D„) and det(G'~^) are positive for m > 0. In more detail, I will 
prove that, like det{Dstag), also det(D„) is real and strictly positive for m > 0. It follows 
from Eq. 02.121) that det(G~^) is real positive (for the issue of zero eigenvalues of G~^, see 
Sec. HV^ . 

I begin by noting [l5| that [75 ® Csj-Dnbs ® ^5] = D^.^^ It follows that [75 (g) ^sj-Dn is 
hermitian, and det{Dn) is real. Moreover, complex eigenvalues of Dn must occur in pairs 
with conjugate values. Therefore det{Dn) will be (strictly) positive if all the real eigenvalues 
are (strictly) positive. 

By "undoing" the all the blocking steps one can express the blocked propagator as 

D-i = i?„ + g„D;4Q|,, (A6) 

where Qn = Q^"^Q^"'~^^ ■ ■ ■Q^^\ and Rn > is determined iteratively from _Ro = and 
Rk = -|- {{z^'^^y /16)Rk-i (for the free theory, see Eq. flH2p ). It is now straightforward 
to show that, if \l/ is an eigenstate of Dn with real eigenvalue A, then 

\=Rn + ^'Qn . ■ (AT) 

^ ^msls '"' 

I have used that Dgtag = D^sis + where the massless staggered operator Dmsis is anti- 
hermitian. For m > 0, it follows that Rn < A~^ < 00. Hence all the real eigenvalues of Dn 
are (finite and) strictly positive for m > 0. 



APPENDIX B: MORE DETAILS ON THE PULL-BACK MAPPING 



Here I discuss in more detail the pull-back mapping introduced in Sec. Ill A[ First, consid- 
ering the original as well as all the blocked gauge fields as a fixed background, let us discuss 
the pull-back mapping Tp'^^ of the fermions only. In analogy with Eq. fl2.14p it is defined 
for — 1 < j < n — 1 by 



k=j+l 



exp 



k=j+i 



(k) 



(Bl) 



The transformation (|A1|) implies — ,^ 
If Dn has no real eigenvalues, the strict positivity of det(£'ri) follows trivially. 



7J if the Dirac and the taste matrices both act from the left. 
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As an example, consider the action of T^"" ^'"^ on a fermion bilinear. Using Eq. (IBip one 
has 



(n— l,n) 



+ 



QW(^(n))^("-i)l L(" '^Q(")^(y("))l . (B2) 



The resemblance to Eq. fl2.10bp is evident. Notice that for x'-"-' 7^ y*^"-* there is no contact 
term. This generalizes to the product of any number of fermion and anti-fermion fields. 
Thus, the fermion pull-back mapping realizes the operator relation 



(B3a) 
(B3b) 



where the right-hand sides were defined in App. [XJ and where the ~ sign means equality 
up to the contact terms that arise when a fermion and an anti-fermion reside on the same 
site of the coarse lattice. Observe that, for the fermion kernels of App. [Aj if no fermion and 
anti-fermion reside on the same site of the coarse lattice, then no contact terms will arise 
under the pull-back Tp'^'' for any j. 

Next let us consider the action of the pull-back mapping on the gauge fields as well. First, 
a few more details on the gauge-field blocking kernels are needed. The non-linear blocking 
kernel Bh is constructed as a sum over the links of the k^^ lattice, 



V 



(k) 



(fe)' 



w 



(B4) 



A simple choice, consistent with the gauge-transformation properties of the fermion kernels 
(1A5D. is 



w 



ik) 



IJ,,X 



(fc) 



(fc-1) ^(fc-1) 



(B5a) 
(B5b) 



where > is a new blocking parameter. Since W G SU{3), one can use the invariance 
of the Haar measure to show that A/fcCl^*-^^^-*) in Eq. (l2.4bD reduces to a numerical constant. 
Many other choices of Bk are possible, see e.g. Ref. |3Cll|. 

Considering the defining equation fl2.2al) . the gauge- field blocking kernel Bk may be viewed 
as a generalized action. This generalized action couples each A; -lattice link to the 

gauge field on the {k — 1)*^ lattice, but it does not couple directly any two fc*^-lattice links. 
As an example, let F{g) denote some function of g E SU{3). Considering an operator of the 



form -^(V^^'^ln)) where O*^"^ does not depend on V^Vf^„), it follows that 



in) 



fl,X'- 



/^-(n— l,n) 



q-{n-l,n) p fyin) \ rj-{n-l,n) Q{n 



(B6) 



Returning to the general case, while the explicit expression gets more complicated with 
every pull-back step, the pull-back mapping is ultra-local because the blocking kernels are. 
If O^"^^ has a compact support, then the support of T*^-''"^(9*^"'' will only slightly increase for 
all j < n. 
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APPENDIX C: ENSEMBLES OF BLOCKED CONFIGURATIONS 



Having introduced all the blocking kernels (see App. |A]and App. [B]), let us discuss the 
generation of ensembles of blocked gauge fields. This is necessary, for example, for the com- 
putation of the coarse-lattice observables numerically. The issue is how to generate blocked 
gauge-field configurations from pre-existing fine-lattice configurations. Using Eq. (12.121) in 
Eq. (112D for t = 1, we get 

/n n 
VU W [PV^'^)] exp ( - - 5^ icf) det"'(D,i,,) (Cla) 
fc=0 fc=0 

= j VUexp{-Sg) dee^ iDstag) = Z{nr) . (Clb) 

Equation (IClbp reminds us that, in the original staggered theory, the Boltzmann weight 
of the fine-lattice gauge field U has the form exp{— Sg) det""' {D stag) with rir = |, |, or 1. 
Equation (IClaP provides a Boltzmann weight for all the blocked gauge fields as well. In view 
of Eq. (IClbp . the process begins with an ensemble of fine-lattice configurations generated 
in the usual way. Given a "mother" configuration Ui in this ensemble, one can generate a 
"daughter" configuration of the once-blocked gauge field vj;^^ by a new Monte-Carlo process, 
by taking exp[— /C^"* (V^^-* , Ui)] as a Boltzmann weight while holding the fine-lattice gauge field 
Ui fixed. As it should, the probability to obtain the pair {Ui,Vi^^} is given by the product 
of the original (normalized) Boltzmann weight (rir) exp{— Sg(Ui)) det"^'' (DgtagiUi)) and 

the new Boltzmann weight exp[—JC^^\vj;^\Ui)]. The process may be repeated on further 
blocking steps 1 < k < n, each time generating a /c*'^-lattice daughter configuration from 
the existing {k — l)*'^-lattice daughter configuration vf'^ using the Boltzmann weight 

exp[-/c(J-)(V(^'),Vf-^))]. 

Once an ensemble of fine-lattice configurations and of daughter configurations for all 
< A; < n has been generated, it can be used to calculate any observable. The blocked- 
lattice fermion propagator can be calculated by repeatedly applying Eq. (I2.1UI) until an 
expression involving the fine-lattice propagator D~^\g is obtained. The blocking kernel Q^^'^ is 
an explicit functional of V'-'^"^^ only, therefore it should be evaluated using a (/c — l)*'^-lattice 
daughter configuration.^^ 



APPENDIX D: LATTICE SYMMETRIES UNDER THE BLOCKING TRANS- 
FORMATION 

As explained in Sec. Ill C[ thanks to the pull-back mapping each coarse-lattice observable 
is at the same time a fine-lattice observable; as such, it is constrained by all the staggered- 
fermion symmetries of the original theory. In a more technical sense, a given fine-lattice 



In the limit [3k ^ oo the blocking Boltzmann weig ht exp[-/C)3^^(V('=\v('=-i))] collapses to a 5-function, 
and the daughter configurations reduce to well-defined functionals of the original fine-lattice gauge field. 
In the case of Eq. (jBSp . for example, one finds that that V'^^^-^o) is equal to C/^_2x(o) C^/i,2i(o)+/i: ^^id ^iil*) 
is equal to v''^[r-]\v'"^~]l, , 
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symmetry may or may not survive as a manifest symmetry on the coarse lattice. I will now 
discuss the fine-lattice symmetries one by one. 

Translation and gauge invariance have a prominent role, and they are secured by construc- 
tion. With the fermion and gauge-field blocking kernels introduced above, the blocked-lattice 
action S'„ (cf. Eq. f l2.2bp ) retains these symmetries manifestly. 

The situation is more subtle with respect to hypercubic symmetry. As explained in 
App. 1^1 the fermion blocking kernels transform non-trivially under 90° rotations. As a 
result, with the blocking transformations as introduced in Eq. (12.21) . in fact Sn is not invariant 
under hypercubic rotations. On a closer look, the reason can also be understood as follows. 
Consider the pull-back r(-^'")C(") of some operator from the coarse to the fine lattice. Under 
a fine-lattice rotation, T^~^''^^0^^^ transforms in the usual way. But, because of the non- 
trivial transformation properties of the blocking kernels, the rotated fine-lattice operator 
cannot be obtained as the pull-back of any coarse-lattice operator! In other words, the 
observables of the coarse-lattice theory do not constitute complete representations of the 
hypercubic group. 

This difficulty can be solved by allowing the blocking kernels to depend on additional 
degrees of freedom, or disorder fields. Each blocking step in Eq. (12.21) is promoted to a 
coherent superposition of block transformations summed over all values of the disorder fields. 
The details are given in App. [El Briefiy, a disorder field allows for parallel transporting of 
the fermion variables of a given 2^ hypercube to any of its sixteen sites in turn. Another 
disorder field allows for all possible orderings for traversing the axes. The gauge-field blocking 
kernels are similarly adapted. With the disorder fields in place, the blocked action becomes 
manifestly invariant under hypercubic rotations. 

Next, the U{1)^ symmetry of the (massless) one-component formalism turns into a 



Ginsparg-Wilson-Liischer (GWL) chiral symmetry [26|, |43| in the blocked-lattice theories. 
See Refs. [H, M for a detailed discussion of both the massless and the massive cases. For 
some further observations, see App. [El 

The last symmetry of the (one-fiavor) staggered theory is shift symmetry. It is generated 
by four anti-commuting elements L^. The action of involves a one- unit translation in 
the /i direction, as well as the multiplication of xi^) and x(^) by sign factors. In the low 
energy limit, shift symmetry reduces to a discrete subgroup of taste-t/(4). The importance 
of shift symmetry is that, without it, a cutoff-scale mass term that breaks the U{A) taste 
symmetry may be induced. This unacceptable mass term is indeed generated if one couples 
the free taste-basis Dirac operator directly to a gauge field [5J], as was found by an explicit 
one-loop calculation in Ref. 53| . 



In contrast, the taste-basis representation constructed in this paper avoids this problem. 
To see this, consider the expectation value of the pulled-back fermion propagator (compare 
Eq. ([M]) )^^ 

G(x("), = ^T(-i'")[V^"^(a;^"^) ^. (Dl) 

Here, cf. Eq. (I2.15p . the subscript "—1" refers to the expectation value of the pulled-back 
operator in the original staggered theory. It is straightforward to verify that the presence of 



Obviously, the size of the translation group gets smaller with each blocking step. 

We may either consider the expectation value in Eq. (|Dip in a fixed gauge, or replace it by a gauge 
invariant one, obtained e.g. by connecting the coarse-lattice fermion and anti-fermion by a coarse-lattice 
Wilson line. 
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the offensive mass term in the coarse-lattice propagator would imply that G{x^"\y^^'') is not 
invariant under shift symmetry [sl, [32|. This is impossible, however, because G{x^'^\y^^^) is 
a correlation function of the original staggered theory. 

As a last exercise that nicely exhibits how the coarse-lattice observables are constrained 
by the fine-lattice symmetries, let us examine the two-point function of a coarse-lattice 
operator with the quantum numbers of an exactly massless (taste non-singlet) Goldstone 
pion [20j], in a tw7o-fiavor theory. As an interpolating field we may take 

vri?(£("^) = V^l'^^x^")) [75 ® ^5] , (D2a) 

where again 75 and ^5 act on the Dirac and the taste indices respectively, and a,b = 1,2, 
label the staggered flavor. Then 

Wr^(n))^W(^wA (D2b) 



1Ti2 (X 



= tr ( [75 ® 6]^n'(5^"\y^"^)[75 ®e5]/^;'(y^"U^"M. (D2c) 



where 



The expectation value on the left-hand side of Eq. (]D2bp is with respect to the partition 
function in Eq. (12. 2p . while on the right-hand side it is with respect to Eq. (IClap .'^^ Via the 
pull-back mapping, the interpolating coarse-lattice fields we use represent specific smeared 
sources on the fine lattice. Now, one can show that Q{x^"'\y^"'^) is strictly positive [isl, 20 



This rules out the possibility of destructive interference caused by these smeared sources; 
the asymptotic decay rate of the correlator must be dictated by the lightest excitation of the 
original staggered theory in that channel, the Goldstone pion. Once again, this shows that 
no fermion mass terms that contradict any of the symmetries of the original staggered theory 
could be generated, because such a mass term would completely change the long-distance 
behavior of this correlator. 

A feature that may be confusing on first acquaintance is that the limiting n —>■ 00 coarse- 
lattice theory is invariant only under 90° rotations, and not under continuous rotations. 
This can be understood as follows. Via the pull-back mapping, even a nominally scalar 
(or pseudoscalar) operator on the coarse lattice has in effect some internal structure for its 
support, that "remembers" the orientations of the axes of the lattice. The lack of manifest 
invariance under continuous rotations in the coarse-lattice theory is, once again, because 
its observables do not constitute complete representations of this symmetry group. In a 
formal sense, the continuum-limit observables accessible by the coarse-lattice theory form a 
discrete subset of the set of "all" continuum-limit observables. If we would keep decreasing 
the coarse-lattice spacing, we expect that the breaking of continuous rotational invariance 
should go to zero like some positive power of QcA. 

In more detail, consider as an example the n — 00 limit of the pion two-point function 
in Eq. flD2b|) . 

^(x,?/;aj = lim (vria^x) 4?(y)) • (D3) 



•^^ In other words, the right-hand side is to be evaluated on an ensemble of blocked configurations, cf. App. [Cl 
•^^ In the free theory, one can check that the operator Drg obtained in the n ^ 00 limit (Eq. (|2.18p ) indeed 
has only hypercubic rotation invariance 2C | . 
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At large (cuclidcan) distances, where the correlator is dominated by the Goldstone-pion 
intermediate state, one expects the factorization 

g{x, y; Oe) ^ e-"^-l^-«l JP-'(n^, aem,) . (D4) 

Here \x — y\ is the usual euclidean distance, and is the unit vector pointing in the 
direction oi x — y. Power corrections that depend on \x — y\ have been suppressed. The 
(direction dependent!) form factor J-{nn, acTn.,^) accounts for the coupling of our coarse- 
lattice interpolating field to the pion intermediate state it creates. The smeared fine-lattice 
operator that corresponds (via the pull-back mapping) to the coarse-lattice interpolating field 
is manifestly invariant under 90° rotations only. As noted above, in the limit Ocm-^ — > 0, the 
form factor T{n^, cLc'^'k) should approach an (n^-independent) constant. 



APPENDIX E: HYPERCUBIC SYMMETRY AND DISORDER FIELDS 

Here I discuss how manifest hypercubic invariance is recovered by summing over disorder 
fields at each blocking step. I will discuss mainly the A; = step, which produces the 
transition from the one-component basis to the taste basis. Subsequent blocking steps work 
essentially in the same way, except that they are somewhat simpler because the subsequent 
blocking kernels act trivially on the Dirac and taste indices. Within this Appendix, I will 
usually use the term "coarse lattice" for the taste-basis lattice obtained via the A; = 
blocking step, in which case I will drop the corresponding superscript-label of the fields and 
the coordinates. 
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FIG. 1: Rotations. The two-dimensional example shows how fine- and coarse-lattice rotations are 
related. The small circle marks the origin. The point X = {\,\) is marked by a cross. Thick 
squares show the blocking pattern. Left panel: Counter-clockwise 90° rotation about the point X. 
The blocked squares are mapped onto themselves: centers are mapped to centers; corners undergo a 
rotation with respect to the square's center. Right panel: The same effect is achieved by a rotation 
about the origin, followed by a translation that brings the point X back into its original position. 



The guiding principle is that we want to preserve the embedding of the coarse (taste-basis) 
lattice into the original fine lattice, i.e. we want to maintain the same breakup of the fine 
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lattice into 2'^ hypercubes. On the coarse lattice, the 90° rotation will be around the origin. 
As can be seen from Fig. [H this corresponds to a fine-lattice rotation around the origin which 
is either followed by, or preceded with, a one- unit translation. Recall that, for staggered 
fermions, the fine-lattice translation group is generated by the four anti-commuting shifts 
that involve a one-unit translation and the multiplication of the staggered fields by sign 
factors. Let us denote by TZ the matrix that produces a 90*^ rotation around the origin. This 
defines a linear, homogeneous mapping of the coordinates, as well as of four-vectors. The 
coarse- and fine-lattice coordinates rotations are given by 

x — ^ x' = Tlx , (Ela) 
x-^x' = n{x) = 7^x + A , (Elb) 

where A is the one-unit translation that follows the rotation (Fig. [1]). The inverse of 
Eq. flETbj) is 

7^-^(x') = n-^x' - A) = n-^x' - n-^A . (E2) 

Making the vector index explicit, the fine-lattice rotation is given by x'^ = 7l{aT)fiu + 5o-^. 
The rotation matrix 7l(aT)^J.u = Sau^riJ. — ^afi^ru + P(aT)nv produces the "counter-clockwise" 
rotation in the (a, r) plane, whereby x'^ = x^ and x'^ = —Xr- Here P[aT)^iv = — ^afi^au — 
Sr^Sri, is the projector on the d — 2 invariant coordinates. For the same rotation one has 
A = 6" in Eq. flElbp . namely, the follow-up translation is in the positive a direction. 

Space-time transformations act on fields by prescribing their value at a point in terms 
of their value at the source of that point (under the "active" coordinates transformation). 
For a coarse-lattice rotation around the origin, the value of the transformed staggered field 
at a fine- lattice point x' will be determined in terms of its value at TZ^^{x'). Therefore, the 
transformation applied to the staggered field is first a shift from x' back to x' — A, and 
then a fine-lattice rotation (around the origin) back to the original orientation, cf. Eq. (lE2p . 
Performing this combined transformation using the rules given in Ref. [8] and plugging the 
result into the right-hand side of Eq. (lAip yields the taste-basis transformation rule for 



hypercubic rotations [53 



V^(x) ^ V^'(x) = [i^®T']V^(7^-'x). (E3) 

Here -R(o-r) = 2~^/^(l — 7o-7r) is the usual Dirac rotation, while T(^aT) = 2~"^/^(7o- — 7r) = ^(L-t) 
produces the rotation on the taste index.^^ 

The fermion blocking kernels flA2p and flASP are already gauge covariant. 1 now discuss 
how to "covariantize" their transformation properties under hypercubic rotations. The idea 
is to simply let any element of the fine lattice, be it a site or a link, transform as it should 
under the above fine-lattice rotation. For the fermion blocking kernels we need to make 
two choices. What choice is being made will be prescribed by a set of discrete-valued 
"disorder fields," that reside on suitable elements of the coarse lattice. In detail, we have to 
decide to which one of the hypercube's sixteen sites will all the fermion variables be parallel 
transported. The chosen site will be determined by a vector field p^. The possible values of 
p^{x) are zero or one, and the fine-lattice coordinates of the chosen site will be 2x + p{x). 
We must also decide in which order to traverse the axes. The ordering will be determined 



■^^ The transformation rule in Ref. [53] fooks slightly different due to a further change of basis. See also 
footnote EH 
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by another coarse-lattice field Cj = lj{x) that takes values in 5*4, the permutation group of 
four elements. 

Let us next establish the transformation rules of these fields. For the chosen-site field we 
demand that, if x = 2x + p{x), then this relation will be respected by the rotation. With x' 
and x' given by (lEip . we must require x' = 2x' + p'(x'). It is straightforward to show that 
the required transformation rule is 

p'(x') = 7^ p{n-^x') + A . (E4) 

Let us verify that Eq. flE4l) is a consistent transformation on this field. We must verify 
that for = 0, 1, also p'^ takes only these two values. Let us again consider the counter- 
clockwise rotation in the {a, r) plane. Only the a and r components undergo a non-trivial 
transformation, which reads exphcitly p'^{lZ(^cjr)x) = Pa{x) and p'^{1Z(^ctt)x) = —pr{x) + l. We 
see that the translation by the unit vector +cr acts precisely to bring p^(7?.(crT)5) back into 
the allowed range. 

In order to write down the transformation rule for the axes-ordering field, let us use 
the defining representation of the permutation group 5*4 in terms of four-by-four orthogonal 
matrices, each of which has one entry equal to one and the rest equal to zero on every 
raw or column. The axes ordering is then given by letting this matrix act on the constant 
four- vector v = (1234), that is, act on the four- vector whose entries are given by = p. 
With this, the transformation rule is 

uj'{x) = Cj{TZ-^x') 7r(7^) , (E5) 

where the permutations %:(JZ{„r)) ^ 'S'4 is represented by the four-by-four matrix 7r(7^(o-r))/^j/ = 

For completeness, recall the transformation rule of the fine-lattice gauge field, which is 
conveniently expressed as [8| 

U{x, y) U'{x, y) = U{n-'x, n^'y) , (E6) 

where 

{u^j,,x , y = x + p, 
Uly,y = x-p, (E7) 
, otherwise . 

We are now ready to introduce new parallel transporters that transform covariantly under 
rotations. Let 

W{x,y,7r;U), (E8) 

be the parallel transporter from y back to x, which traverses the axes in the order determined 
by vr G S'4, as follows. With the constant four vector v introduced above, we let f^- = vrw . 
Starting at y and letting u = ^'^^(4), we first go along the axis until the coordinate 
is equal to x^. The direction is determined by the sign of yi, — x^. Then we go along the 
axis specified by u^(3) and so on. In four steps, each involving a straight line, we go from y 
back to X. Note that the parallel transporter in Eq. (lASp corresponds to the special case of 
choosing vr as the identity element. 

Armed with the more general parallel transporter flE8l) . we modify the fermion kernel 
of the k = step by replacing W of Eq. flA3|) with >V(2x + p{x),2x + r{x),u{x)]U). For 
convenience, the dummy summation variable of Eq. (IA2p has been promoted to a field; its 
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transformation properties are, obviously, the same as those of p{x). The so-constructed 
parallel transporter transforms as 



W{2x + p{x), 2x + r{x),Lj{x)]U) 

= w{2n-^x + piTz^^x), 2n-^x + r{n-^i),Cj{n-^})-u) . 



(E9) 



With this, the right-hand side of Eq. flA2p attains the same hypercubic transformation 
properties as the taste-basis field, cf. Eq. (]E3p . 

Now that the fermions are each time parallel transported to a different hypercube's site, 
we must also modify the gauge-field blocking kernels, so as to maintain gauge invariance of 
the coarse-lattice theory. For the /c = step, the new gauge-field blocking kernel is obtained 
by replacing Eq. (]B5bp with 



fJ,,X 



W{2x + p{x), 2{x + fi)+ p{x + /i), W) 



(ElO) 



where again W is defined by Eq. (lESp . In Eq. (lElOp . the axes ordering is chosen indepen- 
dently for each coarse-lattice link, according to a new disorder field fi^(x) taking values in 
the permutation group S^. Introducing notation analogous to Eq. (lE7p . 



{fi^(x), y = x + p, 
, otherwise , 
its transformation rule is (compare Eq. flE5l) ) 

n{x, y) n'{x, y) = n{n-^x, n-^y) 7r(7^) . 



(Ell) 



(E12) 



We are now ready for the implementation. A complete set of disorder fields is introduced 
at each blocking step. Re-instating the blocking-step label, the "measure" for the disorder 
fields is 



(k) 



24 ^ 



(Ii('=)(x(fe))e54 



n 



(fe) 



16 ■ 24 



E E 

p'-J^\xW)=0,l hi^\xW)eSi 



(E13) 



The k = blocking step takes the form 



VUVxVx eM-SM-xDstag{U)x] 



X 5^ / exp [ - So (V(°),W) -Mo{U) 

(0) 

. r{0) 

X exp 



(E14a) 
(E14b) 



X?V(o)I)7/;(°)DV^^°^ exp 



(E14c) 



The blocking transformation is introduced on line flE14bl) . As promised, it consists of a 
coherent superposition of "elementary" blocking transformations, each corresponding to a 
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particular set of values of all the disorder fields. In going from Eq. (lE14bp to Eq. (lE14cl) we 
both integrate out the original fields, and sum over all values of the disorder fields. Similar 
coherent superpositions of blocking transformations are introduced in all subsequent steps. 

With all the disorder fields in place, the coarse-lattice theory obtained at the n*'^ step 
is manifestly invariant under hypercubic rotations. To see this, observe that all the block- 
ing kernels in Eq. fl2.3p become hypercubic-rotation invariant thanks to the transformation 
properties endowed to the disorder fields. The original action is invariant too, and the sum of 
the original action plus the blocking kernels may be regarded as a generalized action, which 
is hypercubic-rotation invariant as well. The effective action Sn obtained after integrating 
out any number of fields retains the same invariance. 

Because the values of the disorder fields can vary locally, the resulting action Sn can be 
written as a sum over coarse-lattice sites of a local hypercubic scalar. The same would not 
be true had we restricted the disorder fields to take globally constant values only. Such a 
global sum would amount to averaging correlation functions of different (blocked) theories, 
and the result would in general violate clustering. If (manifest) hypercube symmetry was 
enforced by global averaging, violation of clustering would occur at every blocking level k, 
where it is expected to scale like a power of (and, ultimately, like a power of Oc), which 
is unacceptable. This unpleasant situation is avoided, however, because the disorder fields 
are local fields. 

Turning to the representation (12.81) . in order to maintain the gaussian nature of the re- 
maining fermion integral one must refrain from integrating out any fields other than fermions. 
This means that summations over the disorder fields must not be carried out explicitly as 
well. The representation then takes the form 



/I L 
w m k=o 

X J exp - , (E15) 

Related subsequent equations (e.g. in Sec. Ill El) are modified accordingly. We should also 
modify the process of generating ensembles of blocked-lattice configurations. Integrating out 
the remaining fermion fields as well as all the blocked gauge fields in Eq. (lElSP we obtain 



Z = Y.Y.---Y. [ ^^Wexp(-S',)det(D,i,,), (E16) 

(0) (1) (n) 

which is to be compared with Eq. flClbl) .'^^ This equation states the (obvious) result that, 
with no more blocked gauge fields around, the disorder fields decouple from the original 
theory. Therefore, the original fine-lattice gauge field is to be generated as always with 
its usual Boltzmann weight exp{—Sg) det{Dstag), while all the disorder fields are to be 
generated with a fiat measure. Any "tensor-product" configuration made of a fine-lattice 
gauge-field configuration and a configuration of all the disorder fields then serves as a mother 
configuration for the production of the chain of daughter configurations of the blocked gauge 
fields. In practice, this entails the simple instruction that a new set of values of the disorder 
fields is to be picked up at random for any new evaluation of a blocking kernel. 
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The generalization to 7^ 1, cf. Eq. (jCip . is straightforward 
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Last let me address the following question. The fermion blocking kernels are not invariant 
under ordinary chiral transformations, and this leads directly to the GW relation, and to 
the replacement of any ordinary chiral symmetry by its GWL cousin, as discussed for the 



case at hand in Refs. |l5l. |20[|. In comparison to hypercubic rotations, the invariance (of the 
massless limit) under ordinary chiral symmetries is in fact lost already in the free theory 
after one or more blocking steps. One may wonder whether a similar trick with some new 
disorder fields would help us retain the invariance under ordinary chiral transformations in 
the blocked theory. The answer is yes, but it carries with it very little gain as I will now 
explain. 

In the free theory, the generator of (ordinary) chiral transformations in the taste basis 
is [75 ^5]. Within the blocking process, we may enforce the invariance under the same 
(global) chiral symmetry by augmenting each of the fermion blocking kernels with a new 
disorder field M^''^ = M'^'^\x^^'>) transforming like a mass spurion. Specifically, M'^^\x^^'^) 
is a sixteen by sixteen matrix, labeled by a double, Dirac and taste, index. It takes values 
in the U{1) group exp(i6'[75 ® ^^]). Again taking the = step as an example, the new 
blocking kernel would take the form 

ICf = «o(V^^°^ - XQ^'^^KU)) M(°) (^W - Q^'\U)x) , (E17) 

where the dependence on all other disorder fields has been suppressed. Under a global U{1) 
chiral rotation, and assuming that the taste-basis field transforms as — > exp(z6'[75 ® ^5]) ip, 
the new disorder field transforms as 

M('^)(£('=)) exp(-2z0[75 6]) M('=)(£(^')) . (E18) 

Taken alone, the new blocking kernel (1E17P is in fact invariant not only under global but 
also under local U{1) chiral transformations. However, the original action is only invariant 
under the corresponding global transformations (the U{l)e symmetry) in the massless limit. 
Hence, the blocked action obtained after integrating over the original staggered fields, as 
well as over all values of the new disorder field, will be invariant under global chiral rotations 
only, as it should. Of course, that blocked action is not bilinear in the fermion fields, nor 
can it be reasonably approximated by any bilinear fermion action even in the free case. This 



has to be so, or else the Nielsen-Ninomiya theorem [55| would be violated. 

More relevant is the role of the chiral disorder fields within the representation flElSp . 
Particularly illuminating is to consider their effect on the fermions pull-back mapping fIBip . 
It is easily seen that the role of M^^\x^''^) is to multiply the contact term in Eq. flB2p . 
obtained while undoing the k^^ blocking step, by a chiral phase which depends on the Dirac 
and the taste indices. Under a "complete" fermion pull-back all the way to the staggered 
theory on the original fine lattice, nothing would depend on the chiral disorder fields, apart 
from the contact terms encountered along the way. The integration at each blocking step 
over all values of M^''\x^''^) for all x*^^^ would thus wipe out all the contact terms. 

For non-coinciding points x*-"-* 7^ y*-"-* on the "last" coarse lattice, however, contact terms 
are absent anyway. The upshot is that, when evaluating blocked fermion propagators on 
blocked ensembles (cf. App. [U]), we have the following choice for coinciding coarse- lattice 
points. We may either evaluate the contact terms generated by the fermions pull-back 
mapping, assuming there were no chiral disorder fields; or else, we are free to drop them, 
assuming that these disorder fields were present. Whichever choice we make, it is of no con- 
sequence for any non-coinciding coarse-lattice points. The coarse-lattice fermion propagator 
between non-coinciding points is independent of the chiral disorder fields. 
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APPENDIX F: NON-LOCALITY OF THE INTERACTING THEORY AT NON- 
ZERO LATTICE SPACING 



In the free theory, a local square-root operator may be constructed at non-zero lattice 
spacing in the massive case [5Q]. This is not possible in the interacting case: the fourth- root 
theory, or the square-root theory for that matter, are non-local for any non-zero fine-lattice 
spacing aj. This paper shows that the magnitude of all the non-local terms, (but not their 
range!) goes to zero with the fine-lattice spacing. 

Here I briefly repeat the argument why, in the interacting fourth-root theory, the range 



of the non-local terms must be a physical scale [15|] . With Eqs. fl2.8p and (12.121) in mind let 

us assume on the contrary that, after n + 1 blocking steps, a local fourth root exists in the 
sense that ^ 

deti/^(D„) = exp(-i5Se//) det{D ) , (Fl) 

where SSeff is local, and where D is a local lattice Dirac operator which describes one quark 
in the continuum limit. '^'^ 

Let us now compare the actual Goldstone-boson spectrum of the ordinary staggered 
theory (no roots) to that dictated by Eq. fIFip . Substituting the fourth power of Eq. fIFip 
back into Eq. (I2.20p and noting that det^(l)) = det{D ® 1), the assumed locality of 6Seff 
would imply that the RG-blocked theory is in fact a local four-taste theory with an exact 
U{4) taste symmetry. This would imply, in turn, that the fifteen pseudo-Goldstone pions 
must be exactly degenerate. This conclusion is wrong, however! As explained above (see in 
particular Sec. Ill Gl and App. [D]), the RG-blocked theory has the same low-energy spectrum 
as the original staggered theory. This spectrum constitutes of fifteen non-degenerate pseudo- 
Goldstone pions for any non-zero fine-lattice spacing [1, [H, 57]. Thus, the different lattice 



symmetries of the staggered theory and of the putative theory defined by the Dirac operator 
D ®1 rule out a local SSeff- As discussed in Sec. Ill El and on, however, the notion of a 
reweighted theory, namely of a taste-symmetric theory that only approximates the staggered 
theory, can be very useful. 



APPENDIX G: THE PROPAGATOR BOUNDS 

In this appendix I collect a few observations on the propagator bounds (13. 8p and (15.20 . 
The first thing to notice is that the precise form of these bounds is not important, so long 
as it is known that the norm of the inverse Dirac operator in question is bounded by some 
non-zero constant in the limit ay ^0. That constant will depend on mr{ac), and may 
depend on Oc and A as well. 

A configuration for which the bound (13.81) is nearly, but, not quite, saturated is an 
instanton with size p ~ Oc. The bound is not fully saturated because there is no index 
theorem for staggered fermions. (For related observations, see Refs. Ill, [11; for related 



numerical work, see Ref. 29j.) For zero modes of larger-size instantons, the bound (13.80 may 



be corrected by factors of log(p/ac) due to both wave-function and mass renormalizations 



Both SSeff and D may in general depend on the original as well as on all the blocked gauge fields. The 
arguments simplify a bit if no blocking steps are done, as was assumed in Ref. [l5!|. 
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over the range from Oc to p. Because ac < p < A~^, and both Oc and A are held fixed, I have 
neglected such logarithmic corrections. 

My remaining comments concern the bound (15.21) in the reweighted theory. (The bound 
(15. 2p pertains to the one-taste reweighted theory derived from the fourth-root staggered 
theory, but similar comments apply to the four-taste reweighted theory derived from the or- 
dinary staggered theory.) The U{1)^ symmetry of the massless staggered theory is disguised 
as a GWL chiral symmetry by the RG blocking. This chiral symmetry is broken only by the 
(staggered) mass term, which, in turn, is protected against additive mass renormalization. 
The same is not true for the reweighted theories: for fixed aj > 0, the taste-invariant opera- 
tor Dinv,n has no chiral symmetry in the limit where the staggered mass goes to zero. Thus, 
Eq. (15. 3p refiects the presence of an additive fermion-mass renormalization in the reweighted 
theory. 

One should, however, be careful with the interpretation of Eq. (15.31) . First, as noted 
above, already in the staggered theory itself the bound (13.81) is not expected to be completely 
saturated. Second, in this paper I do not consider the reweighted theories as "stand-alone" 
coarse-lattice theories. Their renormalization is defined with references to the underlying 
fine-lattice cutoff aj (see Sec. IIV A2p . Thus, the rightmost term in Eq. (15. 3p vanishes in the 
continuum limit a/ with, in particular, Oc fixed, which implies that no fine-tuning of 
the fermion mass in the reweighted theory is needed. 

The conclusion that an additive mass renormalization of a certain size is present cannot, 
in any case, be drawn based on the magnitude of A„ alone, as can be seen from the following 



example. In Ref. [15| another family of reweighted theories was constructed with a taste- 



invariant Dirac operator Dov,n{f^) = Dov,n{f^) ® 1 such that once again Dn{m) — Dov,n{f^) 
scales in essentially the same way as does Ari(m) = Dnim) — Dinv,nijn)- Nevertheless, 
Dov,n{'n^) satisfies a GW relation in the limit where the staggered mass m goes to zero. This 
imphes that, just like the original staggered mass, the fermion mass residing in Dov,n{jn) 
renormalizes multiplicatively.^^ 

In the physical one-fiavor theory the only chiral symmetry is anomalous, and instantons 
modify the quark's mass. Correspondingly, there is a related tiny correction to the de- 
nominator in Eq. (15.21) . coming from the integration over instantons with size in the range 

< p < etc- By choosing small enough, instantons with size p < ttc are strongly sup- 
pressed, and this correction can be made arbitrarily small relative to mr{ac)- No similar 
correction exists in a theory with more than one degenerate fiavor.^^ For related observa- 
tions on the fourth-root reg ularization of one-fiavor QCD (prompted by the claims made in 



Ref. [491 ), see Refs. lUl, [l 



APPENDIX H: SCALING AND THE MULTI-GAUGE-FIELD REPRESENTA- 
TION 

In this appendix I further expand on scaling issues within the multi-gauge-field diagram- 
matic expansion, discussed in Sec. IIV CI of the main text. This appendix is organized around 
a few examples that each illuminate some particular aspect. In App. IH II I discuss the free 



This observation is relevant for Ref. 



58| where a comparison of the staggered ensemble to a reweighted 



overlap ensemble was attempted. See also Ref. [15|. 
I thank Mike Creutz for a discussion of this point. 
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theory, and in App. IH2l the interacting theory. 



1. Free theory: small-momentum expansion of D„ 



Extending the result obtained in Ref. [20] to m 7^ and adapting to the present conven- 



tions and notation, the blocked free propagator takes the form 

= Mn-Y. {^b^ ® + ® U,]Bl) , (HI) 

where, in the massless limit, 

n 

Mn =i?„ = V(16)*=-7afc. (H2) 

m=0 

k=0 

A straightforward calculation gives^^ 

= ^X^(l + 0(P^)), (H3a) 
B; = 2--1 a, (1 + Oip')) . (H3b) 

Inverting Eq. flHip we find'^^ 

D^{p) = m + i[^®l]+afJ2b5^^5Qpl-Rn(m + i\^®l]y + --- . (H4) 

The ellipsis stand for terms of homogeneity degree three or higher in p^ and m. Observe that, 
even though D„ lives on the coarse lattice, the first three terms on the right-hand side of 
Eq. (lH4p are exactly the same as in the usual taste-basis Dirac operator [sl, IgJ. In particular, 
the "skewed Wilson term" (that comprises the leading taste-breaking term) has a coefficient 
that scales with the fine lattice spacing a/ = 2~"'~^ Qc- Because —n/ac < Pfi < T^/dc, the 
leading taste-breaking term is indeed of order CLf/al. Considering the massless limit we see 
that there is also a term RnP^ with the same structure as an ordinary Wilson term. The 
coefficient i?„ scales with the coarse-lattice spacing because a^^ = 0{ak) by assumption. 
The i?„-dependent terms reflect the fact that Dn satisfles a GW relation in the massless 
limit. 



^•^ Equation (jHSp follows from Eq. (11) of Ref. 20|. The free propagator £)„ ^ is constructed as a sum over 
terms with fine-lattice momentum p + {2tt / ac)k^"\ where the coarse-lattice momentum p is fixed, and 

(n) . (n) 

kii takes integer values such that the fc^ -summation samples all of the Brillouin zone of the original 
fine lattice. Thanks to the suppression provided by the blocking kernels, cf. Eq. (lid) therein, any term 
with 9^ is O(p^) at most. 

Equation pl4)l corrects a mistake in Eq. (3.21) of Ref. [15|. 
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2. Aspects of scaling in the interacting theory 



The scaUng of A„ was derived in Sec. IIV CI for any theory whose partition function can 
be cast in the multi-gauge-field form of Eq. fl4.2p . This includes as special cases the ordinary 
and fourth-root staggered theories, with or without reweighting. Here I illustrate some of 
the "inner working" of the scaling of A„. I first discuss two terms that (should) occur in 
A„, and how their functional form depend on the blocking level n. I then discuss how the 
two terms scale. 

The free-theory result (]H4|) together with gauge invariance requires the presence in A„ 
of a covariant, skewed Wilson term 



O^u^ = ajz-' Yl ^ ® e.^] VW , (H5) 

where, using Eq. (]2.10bp and the overall normalization of the blocking kernels in Eq. (lASp . 
the wave function renormalization factor is 

n 

Zn = l[z^'^. (H6) 

k=l 

The covariant laplacian vjT^ reduces to j?^ + O(p^) in the free-theory limit. The superscript 
notation [n] is meant to remind us that, excepting V*^"\ this covariant laplacian depends 
on the entire "tower" of gauge fields, v|r' = V[r'(W, V*^"-', . . . , V*^""^-*). The explicit (compli- 
cated!) form can in principle be computed using E g. fl2.10p . 



Another example is based on the result of Ref. [32j, where it was shown that the taste- 
basis Dirac operator of the k = step (cf. Eq. fl2.10ap and Sec. IIV C\i contains a term with 
the generic form 

Or=a,^;i°V^J4M»^„„. (H7) 

The notation jFjS = J-'fliU) is a shorthand for (1 — W^yiJA))/ {iafgo), where W^yiU) is a 
Wilson-loop operator (without a color trace), and is the bare coupling; jFjS reduces to F^y 
in the classical continuum limit. In Eq. flH7p all indices except color are explicitly shown, 
and pairs of indices are summed over. The dimensionless tensor Ca/3ijf_iu has 0{go) entries, 
and its explicit form is such that the operator violates taste symmetry. 

The presence of both and Op^ in the taste-basis Dirac operator is dictated by 
shift symmetry, which mixes the leading, dimension-four, taste-invariant part of this Dirac 
operator with taste non-invariant terms of dimension five and higher. The precise form of 
OP and Op^ depends on the covariant "blocking" kernel (lA2p that has been chosen for the 
k = step. While Ref. j32| discusses the taste-basis Dirac operator in the case cto — ^ oo, the 
result flH7p generalizes to ao < cxd. 

How the operator "evolves" with blocking is less certain than in the case of 0^^\ 
where we could appeal to gauge invariance (and the free theory) to determine the overall 



The taste-symmetric propagator usually used in staggered perturbation theory is related to the taste-basis 
propagator by a non-local unitary transformation. 
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normalization. Still, based on the fact that the initial taste-basis operator Dq is known to 
contain Op\ one would expect that Dn contains a similar- looking taste-breaking term, 

C??^ = a/V^i:Vw^S^^^.,., (H8) 

where the coefficients c^apij^iu evolve logarithmically. The operator t]^} has similar prop- 
erties to J^fl, except that now it depends on the "tower" of gauge fields J^]^} = 

^i^i(w,v(°),...,v("-^)). 

Let me now consider the contribution of O^p^ to the taste- violating part of blocked observ- 
ables. The point to make is that any operator of the general form (]H8I) will be suppressed by 
the (effective) gauge- field action in the Boltzmann weight (cf. Eqs. ( 14. 2 p and (14.41) ). The un- 
derlying reason that this works is that the discussion below is nothing but a reconstruction, 
using the terminology of non-perturbative ensemble averages, of the familiar diagrammatic 
argument why taste-violating processes mediated by a hard-gluon exchange are suppressed 
by powers of the lattice cutoff [111, 31. 33l . 

Suppose that a "big part" of J^Jm), which accounts for the gauge-field dependence of 0^p\ 
comes from the fine-lattice gauge field, or from a blocked gauge field V'-'^^ with k <^n. Then 
one can devise gauge-field configurations for which O^p^ is too large. In fact, there exist 
configurations for which (9^^ will be 0(l/a/). However, all such configurations are rare. A 
particular example consists of a fine-lattice vector potential A^^^ with the shape of a wave 
packet whose average momentum is p ~ l/ct/ and whose width is ~ l/oc. Such a vector 
potential is coherent over the coarse-lattice scale; its amplitude. A, can in principle get as 
large as 0(l/aj). Were it to happen, this would give rise to an 0(l/aj) value of O^p' . 
However, large values of A are suppressed. By expanding the fine-lattice gauge-field action 
to quadratic order one finds that the action of the "wave packet" is ~ A^a^/ aj. Therefore the 

average value of A is 0(a//a^). (I have neglected the coupling-constant dependence, together 

with any other corrections that scale logarithmically with the lattice cutoff.) This, in turn, 

implies that O^p^ = 0{af/al) as well. Obviously, this particular "wave packet" configuration 

is further suppressed in the ensemble because it has a limited phase space. But a similar 

conclusion is reached for generic fluctuations of the fine-lattice gauge field, when taking into 

account that these fluctuations are uncorrelated over the coarse-lattice scale. Individual 

(local) fluctuation will typically be 0(l/a/). The random-walk sum of ©((oc/a/)"^) such 

fluctuations is 0(a^/aj), if the sum is over a fine-lattice region with roughly the size of a 

coarse-lattice hypercube. The average value, which is the only thing that a coarse-lattice 

field will be sensitive to, is 0(a//a^). Thus, again, one finds that O^p'' = 0{af/al). Last, in 

the case of the local four-taste staggered theory it is clear that, upon integrating the tower 

of gauge fields, and the fine-lattice gauge field in particular, O^p ^ will induce four-fermi (or 

higher dimensional) terms in the coarse-lattice action Sn, which are suppressed by aj (at 
least). 46 



I now turn to the role of C]^ . Recall that, according to Sec. IIVCI A„ scales like a// 



c 



(up to logarithms) on the ensemble of any theory defined by Eq. (14.21) . However, one cannot 



A parallel statement in the context of the fourth-root theory would only be meaningful within a diagram- 
matic expansion augmented by the replica trick. I stress that the discussion of Sec.|V]is non-perturbative, 
and free of this limitation. 
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deduce that O^-* scales like af/a^ all by itself. We only know that the sum of all contributions 

to A„, coming not only from (9^^ but from many other (higher-dimensional) terms, must 
scale like a//a^. As I will now explain, this scaling depends not only on the suppression 
provided by the gauge-field action (as was the case for O^''), but also on the underlying 
staggered symmetries and, in particular, on shift symmetry. The example below also shows 
that, in other circumstances, there exist operators that would fit the description of O'^^^ as 
given around Eq. flHSp . and yet they would scale as badly as 1/af. 

In order to illustrate in what way things could be different, let us consider the proposal 



5J] to couple the taste-basis fermions directly to a gauge field on the lattice with spacing 



ao = 2a f. The resulting Dirac operator D^^ has no shift symmetry (the superscript ^^DK" 
stands for Dirac-Kahler formulation, which was the main thrust of Ref. 13]). The one- 
loop calculation of Ref. [s^] proves that a taste-violating 0(l/aj) mass term is induced in 
that theory. Thus, (defined in analogy with Eq. fl2.2ip ) will scale like 1/a/ on the 

corresponding ensemble, and will diverge in the limit a/ — >■ 0. In more detail, before doing 
any blocking steps, the 1/a/ divergence originates directly from the skewed Wilson term [53i]. 
By itself, RG blocking obviously cannot "eliminate" any of the divergences of the underlying 
theory. (Assuming on the contrary that a divergent mass term is present in the fine-lattice 
propagator, but absent from the blocked propagator, one reaches a contradiction by invoking 
the pull-back mapping, cf. App. [Dl) Thus, after n blocking steps, the 0(l/aj) scaling is 

expected to come from an operator O^^^^^ with the same general form as in Eq. (]H5|) . Note 
that, because the kinetic term is only marginal, the range of the blocked propagator 
rapidly tends to zero with n. This is, of course, nothing but the decoupling of a fermion 
with a cutoff-scale mass. 

Let us add to the Dirac operator of the DK theory a (taste non-symmetric) mass coun- 
terterm: 

where m^^ = 0{l/af) too; we moreover fine-tune m^^ so that taste symmetry is restored, 
and the correct physical quark masses are obtained in the continuum limit. (Any additional, 
taste-symmetric mass term present in Dq^ will renormalize multiplicatively [a, H^]-) With 
the subtracted operator at the starting point, the taste- violating effects of the DK 

theory have become irrelevant. The desired 0{af/al) scaling of the taste-breaking AfJ^^ 
will now be recovered on the corresponding ensemble. 

The example of the DK theory illustrates that there are two distinct issues here: separa- 
tion of relevant and irrelevant operators; and the scaling of irrelevant operators (in particular, 
within the current RG framework) once we have actually determined what they are. The 
separation of relevant from irrelevant terms will in general require subtractions (i.e. addi- 
tive renormalizations). This is indeed the case in the DK theory. No such subtractions are 
needed in the staggered theory, thanks to its extended symmetry. 

As soon as the appropriate counterterms have been added to the underlying theory, all 
terms in the RG-blocked lattice action that break any of the symmetries of the continuum 
theory must have become irrelevant, and will scale accordingly on the corresponding ensem- 
ble. This amounts to a standard lore in the case of local theories. In this paper I have 



In a strict technical sense, this statement directly applies to internal symmetries. The role of rotation 
symmetry in the RG-blocked action is more involved, see App. |D] and App. [E] 



55 



extended this conclusion to the fourth- root theory (under plausible assumptions). 

In summary, what the example of the DK theory shows is that the only conceivable way for 
the scaling of to go wrong, is when we overlook some of the necessary counterterms of the 
underlying theory. Once all the counterterms needed for the desired (universal) continuum 
limit are introduced, the anticipated scaling of the taste-breaking part of the blocked Dirac 
operator, as discussed in Sec. IIVCI is a generic property of the RG transformation. In 
the staggered case, however, it turns out that no counterterms are necessary jl], thanks in 
particular to shift symmetry. 
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